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Dr  Mazumder  and  Steen  collaborated  at  Imperial  College  during  the  period 
1975-79  when  they  developed  the  first  finite  difference  models  for  heat  transfer 
from  a  laser.  There  were  two  basic  models  developed  that  for  a  stationary 
beam  (ref  1)  and  that  for  a  moving  beam  (ref  2). 

The  predictive  power  of  the  moving  beam  model  was  soon  appreciated  and  found 
to  be  useful  in  helping  with  the  analysis  of  laser  material  processes  such  as 
arc  augmented  laser  welding  (ref  3)  and  laser  welding  (ref  4). 

Mazumder  moved  to  the  States  and  although  both  authors  continued  to  collaborate 
the  distance  between  them  was  a  physical  disadvantage.  This  obstacle  was  partly 
removed  by  the  help  of  this  grant.  As  a  result  the  moving  beam  model  has  been 
greatly  refined  as  described  in  this  report.  The  authors  would  like  toacknowledge 
the  help  of  the  students  in  their  two  research  groups  but  in  particular. 

T.  Chande  (Illinois).  P.  Henry  (ICST)  K.  Lipscombe  (ICST),  M.  Sharp  (ICST). 

They  would  also  like  to  thank  the  USAF  for  enabling  a  useful  cross  fertilisation 
of  these  two  active  research  schools  one  in  USA  one  in  UK. 
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Introduction . 


The  finite  difference  model  on  which  this  work  is  based  is  fully  described  in 
ref  2. 


The  objective  of  this  programme  was  to  develope  this  model  and  use  it.  Thus 
this  report  falls  into  two  parts.  Firstly,  aspects  of  the  development  of  the 


model  to  allow  for: 


1.  Latent  heat  effects. 

2.  Surface  reflectivity  variations 

3.  Mode  structure 

4.  Variable  thermal  properties 

5.  Keyholing 

6.  Fluid  flow  in  melt  zone. 


as  well  as  the  perpetual  problem  of  speeding  up  convergence. 

Secondly  analysis  of  process  performance  il.jstrated  by  the  model  which  can 
not  be  seen  experimentally  because  variables  can  not  be  isolated,  or  in  some 
cases  even  measured. 

2.  Model  Development. 

2.1  Latent  Heat  Effects 

In  general  latent  heat  is  absorbed  in  a  leading  location  of  the  melt  pool  and 
released  in  trialing  location.  Thus  the  overall  heat  balance  is  unaffected  by 
the  presence  or  absence  of  latent  heat.  However,  the  exact  melt  pool  shape  is 
expected  to  be  affected.  Infact  a  longer  melt  pool  should  result  from  allowing 
for  latent  heat.  Any  variation  in  melt  pool  shape  has  a  profound  effect  on  cooling 
rates,  solidification  rates,  and  thermal  gradients.  These  variables,  in  turn, 
describe  the  solidification  mechanism  and  hence  affect  the  weld  metallurgy.  Thus 
to  allow  for  the  latent  heat  is  Important. 

Several  algorlthpts  were  tried  both  at  Illinois  and  London.  The  three  algorithms 
which  were  most  successful  were  based  on  a  variable  specific  heat. 

One  has  an  artiflcally  high  specific  heat  function  during  the  assumed  melting  and 
boiling  temperature  bands  (described  in  appendix  A). 
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The  function  was: 

ce(-r) 

c,(t) 


0.  +  AH  M-  f  •">  T*,  S  T  =£  T„+  at„ 
CPo  +  4Hy/AT,  jtsm  T,  S  T  S  T,  +41; 


where 


the  normal  constant  specific  heat  J/kg  C 


Ah„ah„  = 
at„,at,= 


Transition  energies  -  latent  heat  J/kg 
Temperature  band  over  which  transition  occurs  °C 


although  this  model  was  felt  to  have  good  physical  justification  it  predicted 
lengthy  melt  pools  (  8mm  for  a  2kW  laser  power).  The  recovery  of  all  the  latent 
heat  of  vaporisation  in  the  trailing  edge  of  the  pool  was  considered  to  be  the 
cause.  Improved  versions  of  the  model  which  ignore  the  recovery  of  the  latent 
heat  of  vaporisation  gave  reasonable  melt  pool  lengths. 

Another  latent  heat  algorithm  was  to  have  a  .numerical  hold  for  a  set  number 
of  interations  before  passing  the  transition  temperature.  This  has  an  analogy 
with  the  temperature  rise  of  a  single  point.  For  the  quasi  steady  state  solution 
sort  it  had  little  effect  on  the  answer  compared  to  the  model  without  this 
algorithm 


A  third  algorithm  used  was  to  have  an  artificial  Wnear  variation  of  specific  heat 
with  temperature  to  accommodate  the  latent  heats  over  an  extended  temperature 
range.  For  example:- 


Cr  A  T„ 


CfJT  +  L 


where  A T  =  (T*-T.)  <,«/£=  ««  specific  heat 


used  in  that  interval,  (described  in  appendix  B). 
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This  third  model  gave  the  best  convergence.  The  first  model  had  instabilities 
due  to  the  step  change  in  specific  heat.  The  first  and  third  models  both  had 
resonable  predictive  power. 

This  aspect  of  the  work  is  written  up  in  "Modelling  Laser  Heating  Effects" 
by  P.  Henry  et  al. (appendix  A)and  "Estimating  Effects  of  processing  conditions 
and  variable  properties  on  pool  shape,  cooling  rates,  and  absorption  coefficient 
in  laser  welding"  by  T.  Chande  et  al (appendix  b) 


2.2  Surface  reflectivity  variation. 

The  complexity  in  modelling  laser  interactions  is  the  intimate  interrelationship 
of  the  variables  concerned  .This  not  only  means  experimentation  is  so  complex 
that  a  real  understanding  can  only  be  found  through  mathematical  analysis  but 
also  that  in  building  the  model  it  is  difficult  verify  the  result.  The 
variation  of  surface  reflectivity  with  temperature,  phase  and  surface  shape, 
and  gas  plasma  effects  is  such  an  example.  Hard  experimental  data  at  high 
temperatures  are  not  available  and  if  obtained  of  doubtful  value  due  to  the 
strong  variations  possible  with  such  localised  events  as  surface  deformation 
(including  keyholing)  and  plasma  generation.  Thus  any  algorithm  is  produced 
against  a  blank  background  of  experience.  It  is  known  that  the  material  property 
of  reflectivity  (as  opposed  to  the  surface  property  of  reflectance)  falls  with 
rise  in  temperature  due  to  the  lattice  phonons  and  atomic  vibration  states 
being  more  available  for  absorption  to  electromagnetic  radiation.  However 
the  reflectance  experienced  by  the  beam  is  made  up  of  the  material  property 
reflectivity  , and  a  shape  factor.  The  shape,  in  the  form  of  roughness,  or  melt 
depressions,  may  be  such  as  to  cause  a  single  ray  to  be  reflected  twice  or  more 
before  leaving  the  surface.  In  consequence  greatly  reducing  the  surface  reflectance. 
This  geometric  factor  is  not  temperature  dependent  and  It  may  be  the  dominant 
term; for  example  the  reflectivity  of  aluminium  to  10,6^m  radiation  is  around 
96%  shot  peened  aluminium  has  a  reflectance  of  around  50%. 
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The  gaseous  absorption  in  laser  generated  plasmas  becomes 

significant  at  gas  temperature  of  7000°C  (Saha  Equation  ). 

This  factor  is  highly  temperature  dependent.  The  chosen  algorithm 
needs  to  be  responsive  to  these  main  features  of  the  surface 
reflectivity.  Several  algorithms  were  chosen.  Firstly  a  constant 
reflectivity  up  to  the  boiling  point  after  which  it  was  zero  -  due 
to  the  formation  of  a  blackbody  keyhole.  As  in  the  case  of  latent 
heat  modelling  this  step  function  contributed  to  instabilities. 

Secondly  a  constant  reflectivity  was  assummed  up  to  a  temperature 
T ^ ,  it  then  steadily  fell  to  temperature  T2  after  which  it  was 
zero . 

Using  the  first  model  fig.  10  of  appendix  B  shows  how  insensitive 
the  weld  pool  is  to  the  surface  reflectivity  provided  it  is  low 
enough  to  allow  boiling.  Once  boiling  occurs  the  power  absorbed  in 
the  keyhole  totally  dominates  the  process. 

Using  model  2  with  =  Tg/2  T2  =  Tg  the  same  effect  is  noticed 
where  Tg  =  Boiling  Point. 


2.3  Mode  Structure. 


Considerable  debate  within  laser  manufacturing  circles  is 
centered  on  the  significance  of  mode  structure.  A  Gaussian  beam 
is  known  to  focus  to  the  finer  spot  defined  by  diffraction 
limited  optics  as: 


d 


min 


2 . 4f  X  /D 


where 


f  = 

Focal  length  of  lens 

mm 

X  - 

Wavelength  of  light  used 

mm 

D  = 

Incident  beam  diameter  on  lens 

mm 

d 

mm 

=  Minimum  focal  spot  size 

mm 

In  practice  this  is  never  achieved  because  of  finite  aperture  in 

laser  cavities  and  off  axis  oscillation  states  known  as  TEHmn  mode  I 

(traverse  electromagnetic  modes).  Thus  different  lasers  have 

different  limiting  spot  sizes  for  given  focal  length  lenses  and  j 

spherical  abberration 

I 

d  .  =  2 . 4  (m  +  1)  f X  /D  ' 

min 

I 


where  m  is  a  mode  factor( usual ly  around  4) 

The  model  can  accommodate  any  defined  surface  power  distribution 

expressed  as  the  average  power  over  the  finite  elemental  areas. 

Thus  an  equation  for  the  power  distribution  P  (r,8)  or  P(x,y) 

is  easily  translated  to  the  model  input.  Even  an  experimental 

mode  print  can  also  be  digitalised  for  the  model.  In  trying  to 

describe  high  order  modes  mattiematically  a  problem  is  immediately 

found.  It  is  the  definition  of  beam  diameter.  The  characteristic 

radius  found  in  mathematical  representations  of  various  oscillation 

modes  (eg  ref. 5)  is  a  function  of  the  laser  cavity  alone  and  is 

therefore  independent  of  mode  structure.  What  is  required  for  material 

-2 

processing  is  the  radius  within  which  (1  -0  )  of  the  power  falls  - 

this  factor  makes  the  beam  radius  equal  to  the  radius  for  a  Gaussian 
beam  out  to  the  point  where  the  power  has  fallen  to  /e  of  the 
central  value. 
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s  p  o  1 


R 


C  =  X K 


rdr J©=  (l  -  C 


size  was  defined  as: 

r- 1>  9-o 

Where  P^_  =  total  incident  power  W 

P(r  s)  =  intensity  distribution  W/cmZ 
r,e  =  cylindrical  coordinates 

and  a  mode  factor  K  defined  to  relate  the  processing  radius  R  to 
the  characteristic  radius  of  a  Gaussian  beam,  w 

R  =  Kw 

This  particular  work  is  described  in  appendix  C.  ”An  analysis  of 
the  effect  of  Mode  Structure  on  laser  material  Processing"  by 
M.  Sharp  et  al. 


2.4  Variable  Thermal  Properties. 


The  thermal  properties  describing  the  flow  of  heat  in  a  solid  are 
the  thermal  conductivity,  specific  heat,  and  density.  Together  they 
define  the  thermal  diffusivity.  The  original  model  in  common  with 
most  models  assummed  constant  values  for  these  properties.  However, 
this  variation  with  temperature  in  most  metals  is  by  no  means 
insignificant.  Swifthook  &  Gick  (ref. 5)  used  average  values  over 
the  temperature  range  considered;  but  reference  to  fig.1  below  and 
fig. 1,2, 3  of  appendix  B  shows  this  can  only  be  a  convenience  rather 
thanan  allowable  approximation. 

O 
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define  dun  tn  the  strong  convection  forces  within  r.t.i  molten  pool,  (see 
appendix  1J )  os  well  as  the  lack  of  hard  J.ita. 

Ir  the  model  an  enchanced  thermal  conductivity  of  nrsc  melt  pool  is  assummec 

to  allow  for  convection. 

In  one  version  of  the  model  (see  appendix  B)  the  individual  grid  volumes 
are  given  thermal  properties  by  the  elements  temper,  re. 

According  to  the  algorithm  below: 

K  =  A  +  BT  +  CT2  or  ATn 
Cp  =  ATn 


2 . 1  k e  y  ho  1 i n g 

By  defining  a  grid  element  to  have  boiled  away  if  it  exceeds  the  boiling 
point  4c  allow  the  radiation  to  fall  on  the  element  beneath  it,  a  key¬ 
hole  can  be  computed.  This  simple  model  means  if  keyholing  starts  then 
the  dept  t  u  r  penetration  is  very  large.  Infact  for  certain  power  densities 
it  would  be  infinite  if  the  power  is  sufficient  to  boil  an  element  with 
Vo  r,ui  rsur  d  i  r.p  elements  co  Id.  This  is  not  the  case  in  reality  so  the 
model  had  to  include  some  factor  to  limit  the  penetration  depth.  In 
reality  the  penetration  depth  is  probably  kerbed  by  either  plasma 
absorption  or  keyhole  instability  -  that  is  the  liquid  walls  flowing 
downwards.  (  *pe r imen t a  1 1 y  it  has  so  far  been  impossible  to  seperate  these 
two  possible  mechanism* 

In  the  model  a  plasma  absorption  mechanism  was  assummed  such  that  the 
radiant  intensity  falling  on  the  element  dist  ance  2  from  the  surface 
after  the  elements  above  it  have  toiled  away,  is  defined  by  a  Beer-Lambert 
relationship. 

P  2  -  Pg  exp  (-^  Z ) 


where 


Power  density  at  depth  Z  W/m 

2 


Surface  power  density  w/m 
Beer-Lambert  coefficient 


- 1 
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in-  pr  :■  LI  err-  rmw  shifted  to  tt  t  v.i  [  .•  j!  ii  i  j.ji  ;  jr.ion  wi  Lh  t  ■  :  r  ■ 
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the  values  are  hard  to  use  since  si'.'  r  ■  .  wer-.  net  kn-.wO.  it  ;  : 

possible  to  evaluateyj,  assuming  it  has  a  constant  value,  from 
penetration  data.  Alexander  ref.  6  found  a  slight  plateau  in  weld 
penetration  depth  vs  speed.  Using  this  plateau  penetration  the  following 
values  of p  were  derived. 

(ref  3  Steen)  by  assuming  the  surface  energy  at  the 


Power 

W 

Depth 

mm 

/S 

rn“  1 

Power  at  base 
Keyhole  W 

1  30C 

J-  •  X 

7  00 

260 

1  500 

~i  f 
£.  . 

703 

260 

1700 

2.72 

700 

256 

basis  of  the  keyhole,  whi  reflection,  was  the  same  as  the 

energy  required  to  start  thereyt  , igainst  the  surface  reflectivity. 

The  Beer- Lambert  coefficient  is  a  sensitive  processing  parameter  as 
illustrated  in  appendix  B  fig.  9.  Possibly  the  only  way  a  hard  value 
can  be  found  for  it  is  by  comparing  the  model  with  fusion  profiles. 
This  would  be  a  dangerously  circular  way  of  finding  the  value  of and 
proving  the  model's  validity.  Convergence  is  a  difficulty  when 
keyholing  is  operating.  There  tends  to  be  a  thermal  vibration  on 
certain  grid  elements.  There  does  not  appear  to  be  any  instability  in 
the  heat  transfer  equations. 

The  problem  appears  to  be  associated  with  the  very  high  power  input 
into  the  vaporised  elements  and  the  element  below.  The  vaporised 
elements  are  kept  at  artificially  high  temperatures  in  the  conduction 
matrix.  This  and  the  lack  of  any  sophistication  in  power  distribution 
within  the  keyhole  by  wave  guiding  is  currently  forming  one  area  of 
further  development. 

The  modelling  of  a  keyhole  forming  is  illustrated  in  fig.l  of 
appendix  A. 
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L.b  Fluid  Flow  in  Inc  Melt  Zone . 

The  flow  within  the  melt  pool  would  stimulate  heat  flow  anc  hence  behave 
as  an  augmentation  of  the  thermal  conductivity.  It  differs  from  the 
thermal  conductivity  in  that  it  is  not  a  scalar  property.  However,  in 
the  laser  interaction  zone  the  flow  is  driven  by  a  surface  tension 
force  operating  in  the  direction  of  maximum  thermal  gradient,  and  thus 
on  the  surface  the  flow  vector  is  parallel  to  the  thermal  flow  vector. 
Beneath  the  surface  this  is  unlikely  to  be  the  case. 

The  best  mathematical  model  would  incorporate  a  flow  analysis  via  the 
Nav i er- S tokes  equations.  Such  analysis  would  increase  the  storage 
requirements  for  computing  and  the  computing  time.  However,  a  start 
was  made  on  this  complex  problem  as  discussed  in  appendix  D. 

For  the  thermal  model  an  arbitrary  increase  in  thermal  conductivity 
was  tried.  The  algorithm  was: 

k^  =  m  tk) 

where  m*  2  or  3  for  al 

4  or  5  for  Fe 

The  effect  of  this  modification  was  negligeable.  This  is  probably 
because  an  increase  in  conductivity  of  the  melt  pool  only  should  increase 
the  melt  pool  temperature  and  hence  size  but  the  consequent  increased 
gradient  over  the  outer  area  of  the  pool  undersurface  would  have  the 
opposite  effect. 

2.7  Convergence. 

The  original  model  had  an  exponential  grid  pattern  (ref.1).  This  led  to 
difficulties  in  defining  power  inputs  correctly  and  the  temperatures 
of  elements.  It  was  believed  that  mismatches  of  temp*atures  and  powers 
led  to  mathematical  lack  of  focus  in  the  model  and  hence  some  instability 
Convergence  with  a  uniform  grid  was  comparatively  swift  but  the  boundary 
conditions  were  less  assured. 
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For  example  the  edge  boundary  condition  was  assumed  to  be: 


Thus  the  model  was  arranged  to  cycle  a  fixed  number  of  times  ( 50) 
with  an  exponential  grid  and  then  redefined  to  a  uniform  grid  up  to 
the  limits  of  convergence.  The  removal  from  the  model  of  step 
function  algorithms  as  far  as  possible  also  contributed  to  stability. 

Typical  computing  times  are:-  1500  iterations  in  11  min  55  sec. 

CDC/Cyber  174. 
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3.  Application  of  the  Model . 

The  model  produces  an  answer  which  is  the  temperature  matrix  for  the  quasi 
steady  state  profile  beneath  the  laser  beam.  This  thermal  matrix  is  moving 
through  the  material  at  the  traverse  speed. 

Thus  this  data  needs  adapting  to  whatever  application  is  considered. 

The  usual  parameters  calculated  from  the  model  are. 


a.  Boundaries  of  specific  isotherms  -  eg  fusion  boundary,  or  Ac^  boundary. 

b.  Boundaries  for  regions  above  a  certain  temperature  for  longer  than  a 
certain  time. 

c.  Thermal  gradients 

d.  Solidification  rates 

e.  Cooling  rates. 


3.1  Isotherm  boundary  applications. 

One  of  the  earliest  applications  was  in  the  evaluation  of  the  role  of  the 
electric  arc  in  arc  augmented  laser  welding  by  Eboo  et  al  (ref  3,7)  and 
theoretical  welding  speeds  by  Marzumder  et  al  (ref  4,8)  and  appendix  A  fig 
10,11,12,13.  Currently  appendix  B  shows  the  more  intraverted  applications  of 
assessing  the  significance  of  the  variation  of  parameters  which  are  impossible 
to  experiment  with  such  as  jl  ,  reflectivity,  thermal  conductivity  (fig  9,10, 

11  appendix  B). 


The  experimental  determination  of  depth  of  hardness  was  originally  attempted 
by  Courtney  et  al  (ref  9  )  in  which  the  parameter  P?, 


DV 


where 


P  =  Power  W 

0  =  Incident  spot  size  m 
V  =  Traverse  speed  m/s. 
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was  found  to  be  proportional  to  the  deptkof  hardness  when  there  was  no  surface 
melting.  Analytic  studies  using  the  model  are  discussed  in  appendix  A. 

Fig  14,15  appendix  A  show  that  this  experimentally  observed  linearity  is  only 
approximately  true,  the  beam  diameter  behaves  in  a  more  complex  manner  but  not 
one  causing  sufficient  variations  to  be  noticed  from  the  normal  experimental  noise. 

In  appendix  A  the  weld  penetration  is  calculated  for  large  unfocussed  spots 
because  this  falls  in  the  region  of  rapid  solidification  or  surface  alloying. 

These  laser  processes  are  currently  attracting  considerable  attention  and  the 
model  is  capable  of  reasonably  accurate  melt  pool  size  predictions.  But  it 
must  be  remembered  that  precise  reflectivity  and  beam  sizes  are  not  known  in 
practice. 


3.2  Thermal  gradient  and  cooling  rate  application. 

In  these  processes, as  well  as  welding, the  solidification  structure  is  very 
important.  The  fineness  of  the  solidified  structure  (eg  dendrite  arm  spacing) 
is  found  to  be  Inversely  proportional  to  the  cooling  rate  (G.R) 

(thermal  gradient  G  C/  m  x  solidification  rate  R  m/s)  and 

for  eutectic  growth  "Xr  *  const^ll  ^CSere  "X.  =  plate  spacing, ro 

The  solidification  mechanism  (planar,  cellular,  dendritic)  is  described  by  the 

tendency  to  constitutional  supercooling  which  is  defined  by  the  parameters  (G/R). 

The  phases  which  crystalise  or  undergo  solid  state  transformation  are 

decided  by  thermodynamics  and  kinetics.  During  rapid  quenching  supercooling  is 

normal,  this  can  affect  which  phases  are  favoured  thermodynamically. 

What  actually  occurs  is  determined  by  the  nucleation  and  growth  mechanisms  in  the 
melt.  In  this  respect  the  crystals  at  the  base  of  the  melt  represent  nuclei 
encouraging  epltziial  growth  while  at  the  same  time  undergoing  grain  growth. 

This  extremely  complex  interaction  of  events  In  a  moving  melt  pool  Is  difficult 
to  describe  by  the  data  from  the  model.  The  cooling  rate  G.R  and  some  extra 
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data  . .  j f . f  as  a  ITT  curvt  for  the  alloy  considered  are  currently  the 
ton  f  inalysis  possible,  some  initial  studies  involving  the  basic 
kinetic  and  nucleation  parameters  have  been  made  by  Bergman  (ref/O). 
Whatever  eventually  matures  in  this  area  the  thermal  gradient  G  and 
the  solidification  rate  R  are  almost  certain  to  be  required. 


The  scale  and  extent  of  segregation  in  the  weld  metal  would  be  expected 
to  be  determined  by  the  concentration  "bow  wave"  ahead  of  the 
solidifying  -front  which  if  governed  soley  by  diffusion  in  the  melt 
can  be  described  by 


^ 


f" 


I  -  k 
k 


(-&)) 


where 


k 

X 

D 

R 

x 


partitioning  coefficient 
concentration 
solute  diffusivity  in  the  melt  m^/s 
rate  of  solidification  m/s 
distance  from  solidifying  front  m 


=  X  solid  /  X  liquid 
2 


Thus  a  knowledge  of  R  would  allow  some  semi  quantitative  analysis  of 
segregation . 

Cooling  rates  G.R.  are  discussed  in  appendix  A,B,  and  E 

a 

In  appendix  E  Range  of  rapid  solidification  structures  with  laser 

* 

surface  treatment  of  special  steels  by  G.  Christodoulou ,  P.  Henry, 

W.M.  Steen  Proc  4th  Int.  Conf  on  rapidly  quenched  metals  OlBdai  1981). 


Graphs  are  shown  of  G.R.  and  G/R  contours  in  a  keyhole  weld. 


In  appendix  B  an  advance  is  made  by  identifying  a  dimensionless  melt 
fronf  velocity  0  where  0  * (P/d  2)  /Up  Cp  Aim.) 
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This  is  a  psuedo  melting  efficiency.  This  group  may  prove  to  be  a  useful 
parameter  for  collating  the  data  obtained  from  each  run. 

Dimensionless  operating  parameters  which  allow  generalised  results  to  be 
expressed  is  the  next  step  in  the  application  of  this  model. 

Attempt  to  apply  the  model  to  other  processes  than  welding  are  underway. 
Appendix  F  shows  the  coupling  of  the  thermal  conduction  model  to  a 
variable  surface  shape  and  power  distribution  to  simulate  laser  pneumatic 
claddi ng . 

The  success  of  this  is  quite  striking  (eg  fig  11  appendix  F) 


CONCLUSIONS. 


A  sophisticated  finite  difference  model  has  been  developed  capable  of  taking 
account  of  surface  variations  in  convection,  radiation  and  reflectivity 
bulk  variations  in  thermal  conductivity,  density  and  specific  heat.and 
Physical  events  such  as  phase  transformations  involving  latent  heat 
keyhole  formation,  plasma  absorption  and  flow  in  the  melt  pool. 

The  computing  time  is  reasonable  for  a  large  mainframe  computer. 

Appl icatiorAof  the  model  are  currently  centred  on  calculating  thermal  gradients> 
G^solidification  rates,R,  and  locations  of  specified  isotherms  for  processes 
such  as  surface  melting  and  welding.  It  has  also  been  adapted  to  model  clad 
thickness  and  rate  of  dep  osition  as  well  as  G.R,  and  isotherm  locations  for 
the  process  of  laser  cladding  with  blown  powder  feed. 
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ABSTRACT 

The  results  shown  here  are  some  preliminary  studies  lr,  the  co-ojjer.it  iv.  devr  1  op:  1  :  :  u 

j  finite  difference  model  of  the  laser  heating  process.  Two  models  are  currently  .nb.-i 

investigation  that  derived  at  Imperial  College  and  that  at  the  University  of  Illinois. 

They  have  a  similar  root  in  the  original  work  of  Mazunder  and  Steen  (1)  . 

i.-porte.!  here  are  some  results  of  the  working  of  the  model.  Them  illustrati  a:; 
other  thinns  the  flexibility  of  this  form  of  calculation  and  our  ability  to  disc. on-  •  : 
importance  of  latent  heat,  variable  reflectivity,  variable  thermal  properties  and  key::  line. 

INTRODUCTION 

Laser  irradiation  produces  rapid  local  heating  in  metallic  substrates,  this  can  tie 

followed  by  high  cooling  rates  when  irradiation  is  stopped  or  withdrawn  due  to  conduction 

into  the  substrate.  Metalluraically  attractive  surface  structures  can  thus  h'e  readily 
generated . 

Laser  material  processing  involves  complex  physical  processes  and  interactions  trial  would 
require  considerable  experimental  effort  to  understand.  A  mathematical  model  offers  a 
powerful  addition  to  an  otherwise  purely  experimental  approach  when  it  comes  to  res.  ■  1  vino 
these  oroblens. 

Analytical  mathematical  models  are  limited  by  the  assumptions  necessary  to  make  the 
problem  tractable.  Numerical  models  can  be  more  general  and  can  be  made  to  allow  for  non¬ 
linear  surface  events,  such  as  radiation  losses,  or  variable  absorption  due  to  surface 
evaporation  as  in  "keyholing".  They  can  also  be  made  to  accommodate  variations  in  the 
power  distribution  across  the  laser  beam  (mode  structure)  as  well  as  variations  in  the 
substrate  optical  and  physical  properties  due  to  temperature  dependence. 

A  number  of  numerical  models  of  laser/material  heating  have  been  suggested.  The 
simple  one  dimensional  heat  flow  model  of  Rear  et  al  (2) ,  has  its  attractions  since  the 
answers  can  be  displayed  In  a  fairly  generalised  form,  as  for  the  spot  heating  solution 
of  Ready  ( 3 ) . 

More  sophisticated  three  dimensional  models  are  less  easy  to  express  as  generalised 
solutions.  Kou  et  al  (4)  derive  a  model  based  on  an  oblate  spheroidal  coordinate 
system  which  is  converged  by  a  Gauss  Siedel  iteration  method.  Since  this  coordinate 
pattern  would  be  expected  to  fit  the  isotherms  better  than  the  orthogonal  system 
discussed  here,  it  should  achieve  a  more  accurate  solution  faster.  However,  it  appears 
that  Kou ‘ s  method  involves  a  considerable  computing  penalty,  and  is  not  yet  able  to 
account  for  keyholing  and  surface  heat  losses. 

The  model  described  here  allows  for  surface  heat  losses  by  radiation  or  convection, 
variable  thermal  properties,  and  reflectivity  as  well  as  latent  heat  effects  and 
"keyholing” . 

There  arc  so  many  variables  that  dimensional  analysis  docs  not  significantly  simplify 
the  presentation  of  the  solutions.  Thus  the  working  of  the  models  is  illustrated  here  by 
examples  of  particular  laser/naterlal  interactions. 

DESCRIPTION  OF  MODELS 

In  this  project  two  models  for  laser  material  processing  were  compared.  They  are  both 
based  on  a  3D  numerical  finite  difference  model  described  elsewhere  by  Mazumder  an. I  Mecn  (1)  . 

The  original  root  form  had  an  exponential  grid  to  remove  the  boundaries  so  that  t  in  y 
had  effectively  no  thermal  gradient  across  them.  The  Illinois  model  i  nt  rodu.-.-.l  .  i  w. 
stage  arid  sett  inti,  starting  with  an  cxjioncntiai  grid  to  roughly  find  i  lie  heated  .i.g 
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a,  b,  c  are  constants. 

JVanvhile  the  Ir.perial  College  version  had  been  adapted  by  rboo  (5)  to  ulLow  toi  ...  y- 
holmg  and  plasma  absorption  effects.  They  also  tried  variable  reflectivity  in  lirvu  ;  • - 
f unctions. 

Natural  ly  us  lrprovemen  ts  were  developed  in  «.,•  i  *  r .  *  r  plac;«.  they  vert-  incorporate.:  • 

versions.  Currently  the  analysis  is  focussed  upon  latent  boat  effects.  The  1:. :• 
model  attempts  to  allow  for  this  by  introducing  a  hold  in  the  interaction  of  ^h<-  r-  »>.  t- 
t ion  method  as  transition  temperatures  are  encountered. 


The  Irpenal  model  has  been  adapted  to  have  a  modified,  specific  heat  icnct 
uses  an  artificially  high  specific  heat  during  the  as summed  melting  and  boil; 
The  function  is 

Cp(T)  =  Cp_  *  AH/ AT  fror  T_ <  T  £  T  *  'T 
1  *o  r  n  m  m 

and 

Cp(T!  =  CpQ  *  AHv.'ATv  fron  T,  £  T  *T  (7y  *  ATy) 

where  Cp  =  the  normal  constant  specific  heat  J/kg°C 
o 

AH  , AH  transition  energies,  -  latent  heat  J/kg 
v  r. 


A-m'  ATv 


tenperature  band  over  which  transition  occurs  °C 


The  effect  of  these  differences  and  their  physical  significance  forns  the  nain  area  of 
cooperation  work  between  the  two  centres.  So  far  this  cooperation  has  resulted  in  a  rodcl 
of  increasing  accuracy  and  sophistication  with  decreasing  convergency  tines. 


rig.  1  shows  cross  sections  of  the  nelt 
pool  for  typical  laser  glazing  runs  on  CPTi 
(Commercially  Pure  Titanium! .  The  reflectivity 
in  this  run  changed  linearly  from  a  value  such 
as  C.7  down  to  zero  at  the  melting  point.  The 
onset  of  keyholing  can  just  be  seen  at  the 
slower  speed.  In  this  figure  latent  heat  was 
allowed  for,  similar  shaped  profiles  are  ob¬ 
tained  if  latent  heat  affects  are  ignored, 
however  in  these  cases  the  melt  pool  and 
isotherm  fields  produced  are  wider  and  deeper 
across  the  track  but  shorter  along  it. 

Latent  heat  effects  suppress  the  thermal  pen¬ 
etration  by  introducing  heat  sinks  which  slow 
up  the  heating  process  and  limit  the  peak 
temperatures  and  thermal  gradients. 

In  figs.  2,3,4  the  affects  of  increasing 
traverse  speed  on  the  melting  (1960K!,  1 3CTH 
tOCOK  and  70CK  isotherms,  is  demonstrated. 
Keyholing  is  just  apparent  at  the  slower 
speed  of  18  cm/s,  but  above  this  at  30  or 
ICO  cn/s  the  melt  pool  is  conduction  limited 
with  only  surface  absorption  of  the  laser 
radiation.  The  1CCOK  isotherm  roughly  narks 
the  extent  of  the  KAZ  (heat  affected  zone) 

There  Is  a  noticeably  steeper  gradient  at 
the  faster  speeds.  Figs  5  and  6  are  views 
of  the  melt  pool  fron  above.  The  lines  show 
the  extent  of  the  pool  at  depth  intervals  of 
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Fi'j.  1  Theoretical  melt  pool  cross 
section  in  CI*  Titanium  nivu.ii-  with  l 
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Tru  effect  of  latent  heat 
i tt.  cause  an  onernv  absorption  at  the 
leading  edoe  of  the  pool  and  an  energy 
release  at  the  trailing  edge,  thus 
extending  the  length  of  the  pool.  The 
results  from  the  Illinois  version  which 
does  not  allow  for  latent  heat  by  the 
quasi-specific  heat  method,  but  by  the 
numerical  hold  method  produces  a  shorter 
pool  of  around  1.5  nn  for  these  para¬ 
meters.  It  is  possible  to  argue  that 
the  latent  neat  of  vaporisation  will  be 
totally  or  partially  lost  to  the  system 
and  not  repaid  at  the  trailing  edge  as  is 
assumed  here.  It  is  also  possible  to  aroue 
that  the  reflectivity  is  only  zero  once 
boiling  occurs  and  not  once  melting  occurs, 
both  of  these  arguments  would  produce  a 
shortened  weld  pool  and  variations  such  as 
these  are  currently  under  development. 
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FIG,  1 


SPEED  =30cm/s 


Fig.  2.*4  The  effect  of  increasing  tra¬ 
verse  speed  on  the  700, 'k  1000’k, 
1300'k,  and  1960'k  isotherms  in  CP 
Titanium  using  2  kw  laser  power  and 
0 •  *1  mm  diameter  beam  size. 


Fig.  5.  View  of  melt  pool  from  above 
showing  contours  at  100  urn  depth  inter¬ 
vals:  for  2  kw  laser  power;  0-4  mra 

spot  size  and  30  cm/s  traverse  speed. 


Fig.  6.  As  for  fig.  5  but  at  the  slower 
speed  of  18  cm/s 
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*« - 1 1:  indication  that  there  is  sone 

•  "asria  heating  which  would  create  a 

i  1  heard"  to  the  weld  pool.  The 
: •  : :  ,-n  of  solidification  at  the  trailing 
*e  shows  a  steo»  cooling  rate  of  the 
rdoi  of  10f  -  iO^K/s  in  both  cases, 
with  the  highest  rates  at  the  surface. 

The  "heat  sink"  affect  of  the  latent 
heat  is  illustrated  in  Fig.  9  where 
Fig,  7  is  transposed  from  temperature 
to  temperature  equivalent  based  on 
the  enthalpy 

Teq  =  T  ♦  Iq 

o  s  Any  stored  latent  heat 

It  is  seen  here  that  the  plateau  region 
in  the  temperature  profile  is  actually 
a  region  in  which  the  latent  heat  of 
boiling  is  conducted  out  of  the  system. 

As  argued  previously  this  is  probably  not 
true  physically  and  the  latent  heat  of 
boiling  is  lost  to  the  system  by 
evaporation.  Where  allowance  is  made 
for  this  the  plateau  region  of  Figs.  7 
and  8  would  be  considerably  reduced. 


TIME  [ms  ] 


Fig.  .  Thermal  history  curves  for 
points  at  different  depths  on  the  wold 
pool  cept.ro  line.. 

(Power'.  =  2kw;  spot  size  1  -  0  4mr,l 
(speed  -  100  cm/s  CP  titanium.  ) 


Fig.  8.  As  for  fig.  7  but  at  the 
slower  speed  of  18cm/s 


Fig.  9.  Enthalpy  -  time  curve  of  the 
data  shown  in  fig. 7. 


Figi_14.  A  parametric  plot  of  depth  vs 
P/.DV  for  the  penetration  of  the  1960'k 
melt  isotherm. 

The  symbols  are  for  the  following  con¬ 
ditions  : 


p 

x  -  2kw 

0 

D 

4mm 

V 

various 

♦  -  2kv 

1 

0mm 

" 

C  -  2kw 

2 

Omm 

" 

;  -  Ikw 

3 

0mm 

M 

t  -  various 

0 

4mm 

30cm/s 

Fig.  15  Variation  of  penetration  of 
the  7  00 " k  ,  1000' k,  1100'k,  and  196G'k 
isotherms  with  power 
(D  =  0-4mm,  V  =  30cm/s,  CP  titanium) 


reduced  to  around  0.2  mm.  It  is  thus  possible  to  calculate  the  required  operating 
curves  of  soeed,  power  and  beam  diameter  to  achieve  certain  penetrations  of  specific  iso¬ 
therms  in  any  materials  of  interest. 


Experimentalists  have  tried -to  correlate  penetration  data  against  a  variety 
of  parameter  groups  such  as  specific  energy MfVFV) ,  thermal  penetration  (P//DV)  or  even 
power  density  (P/D2).  Fig.  14  shows  the  theoretical  predictions  of  depth  vs  P//DV. 

A  nesting  of  the  data  is  observed  which  if  it  had  been  experimental  could  be  considered 
as  a  viable  correlation  within  the  bounds  of  experimental  accuracy,  theoretically  we  fir 
a  significant  scatter  to  be  present.  The  fit  to  1//V  is  good  but  that  for  Pvs  1//D 
^"0t.S0  "°°d-  The  reasonfnc1  for  a  Possible  correlation  of  deoth  of  HAZ  with 
P/  DV  is  based  on  the  sinnlif ication  that  the  variation  in  temoerature  with  deoth  can  t 
expressed  in  the  form: 


T*  -  A'  *  where  T*  =  nkD/P  ,  ,, 

and  u*  =  /4ut/x 

The  heating  time,  t,  is  approximately  given  by  the  interaction  time  (U/V)  and  thus 

■i*  =  1  /4uD 

X  V 

substituting  these  terms  into  the  dimensionless  equation  yields 

T"  kD  =  A  /4od 

P  x  — 


and  hence  for  a  given  T  depth  of  penetration  x  is 


/DV  Te"k 


26 


DEPTH  (mm] 


Fig.  1  fa  The  variation  of  cooling  rate 
G.R  with  depth  in  a  keyhole  weld  pool, 
for  the  centreline  and  a  plane  mn 
to  one  side  of  the  centreline. 

(P  =  ikw,  D  -  0-4mm,  V  =  16cra/s) 


Fig.  17  As  in  fig.  H  out  for  a  co 
auction  lir:ited  weld  pool. 

(P  =  2kw,  D  =  0  4mm,  V  =  100  cm. si 


To  fit  the  actual  data  a  more  refined  model  than  equation  (1)  is  beine  developed. 

Fin.  15  shows  the  non-linear  increase  in  depth  of  theoretical  penetration  with  power. 
This  is  partly  due  to  the  temperature  suppression  with  latent  heat  but  principally 
is  due  to  the  rate  of  heat  flow  within  a  solid  body  being  limited  by  the  Fourier  number, 
which  is  not  power  dependent. 

Finally  Figs.  16  (  17  illustrate  the  possibility  of  calculating  important  metallurgical 
oarameters  which  cannot  normally  be  measured.  In  this  case  the  solidification  interface 
cooling  rate  (G.R  where  G  is  solidification  rate  (m/s)  and  R  is  the  thermal  gradient 
k/s)  which  affects  the  scale  of  metallurgical  microstructures  is  shown  for  the  centreline 
and  line  adjacent  to  it  (67  urn  to  one  side).  The  curves  show  the  rise  in  cooling  rate 
as  the  surface  is  approached  as  has  been  predicted  by  Kear  et  al '  and  Kou  and  Mehrabian  '  '. 
Fig.  17  which  represents  a  conduction  limited  track  shows  the  same  kind  of  form  as  that 
predicted  by iKbu.  Kear's  one  dimensional  model  predicts  a  generally  much  steeper 
increase  in  cooling  rate  up  to  an  equilibrium  value.  Fig.  16  which  represents  a  keyhole 
track  shows  an  area  where  the  cooling  rate  is  higher  than  the  material  above  it  this  is 
due  to  the  thfee  dimensional  cooling  of  the  thin  neck  of  the  keyhole.  As  would  be 
expected  higher  cooler  rates  are  obtained  with  the  smaller  melt  volume  of  the 
conduction  limited  track  (an  order  of  magnitude  difference) . 

CONCLUSION 


This  is  a  preliminary  survey  of  some  results  from  the  3D  finite  difference  models  being 
developed  at  University  of  Illinois  and  Imperial  College.  These  models  are  being  developed 
to  minimise  the  number  of  assumptions  which  need  to  be  made  in  the  theoretical  predictions 
of  laser/material  processes. 

Currently  the  models  are  being  used  to  study  the  effect  of  latent  heat  on  weld  pool 
characteristics;  the  correlation  of  operating  parameters  with  effects  within  materials, 
the  mapping  of  important  metallurgical  parameters  (such  as  G.R  and  G/R)  and  analysis  of 
various  physical  models  of  the  laser/material  interaction  of  latent  heat,  variable  reflecti¬ 
vity,  variable  thermal  properties  and  kcyhollhg. 
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APPENDIX  B 


ESTIMATING  EFFECTS  CF  PROCESSING  CONDITIONS  AND 
VARIABLE  PROPERTIES  UPON  POOL  SHAPE,  COOLING  RATES 
AND  A3SORBTIQN  COEFFICIENT  IN  USER  WELDING 

by 

T.  Chande  and  J.  Mazumder  t 
Materials  and  Design  Division 
Department  of  Mechanical  and  Industrial  Engineering 
University  of  Illinois  at  Urbana-Champaign 
1206  W  Green,  Urbana,  IL  61801 

ABSTRACT 

This  paper  examines  the  role  of  traverse  speed,  Beer  Lambert 
absorbtion  coefficient^  ,  surface  reflectivity  and  changing 
liquid  thermal  conductivity  upon  the  shape  of  the  melt  pool, 
and  the  cooling  rates  that  occur.  The  dependence  of  3  upon 
processing  conditions  is  also  examined.  A  three-dimensional 
variable  property,  moving  heat  source,  quasi  steady  state, 
finite  difference  model  for  heat  conduction  into  the  substrate 
during  laser  welding  is  used.  With  an  increase  in  traverse 
speed,  the  pool  flattens  out,  is  swept  back,  and  oooling  rates 
increase.  An  increase  in  fi  sharply  decreases  the  depth  of 
penetration.  With  the  onset  of  melting,  changes  in 
reflectivity  did  not  change  pool  shape  significantly.  An 
increase  in  effective  liquid  metal  thermal  conductivity 
increases  melt  pool  aspect  ratio.  Cooling  rates  increase  as 
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the  energy  density  In  the  pool  decreased.  A  dimensionless 
melt  front  velocity,  0  ,  is  defined  such  that  cooling  rates 
exceed  1000  K/s  as  0  approaches  unity.  The  product  fjz,  where 
z  is  the  depth  of  penetration,  is  shown  to  vary  linearly  with 
the  natural  logarithm  of  0.  These  results  imply  that  0 

i 

affects  depth  of  penetration  more  than  the  width,  that  an 
upper  bound  for  fi  may  be  deduced  from  0,  that  variations  in 
surface  reflectivity  are  less  critical  in  laser  welding,  that 
the  maximum  thickness  that  can  be  welded  in  a  single  pass 
decreases  as  fluid  flow  becomes  more  dominant  in  the  melt 
pool,  and  that  cooling  rates  increase  as  pool  energy  density 
decreases,  especially  for  values  of  0  <  100. 


30 


2.  INTRODUCTION 


Laser  welding  has  received  much  experimental  attention[ 1 ], 
but  not  all  the  physical  phenomena  involved  are  fully 

i 

understood.  Laser  material  interactions  consist  of  complex, 
and  generally  short-lived,  but  intense  events.  Many  important 
aspects  and  effects  of  these  interactions  are  not  directly,  or 
routinely  measureable,  3uch  as  the  absorbtion  coefficient  in 
the  keyhole,  and  the  temperature  and  fluid  flow  velocity 
distributions.  It  is  also  somewhat  difficult  to  isolate  the 
effect  of  a  particular  variable  without  extensive 
experimentation,  due  to  the  many  substrate  properties,  and 
operating  conditions,  contributing  to  the  observed  effects. 

On  the  other  hand,  a  mathematical  model,  based  upon  insights 
gained  during  experimentation,  can  be  used  to  calculate  the 
above  parameters.  A  mathematical  model,  within  the 
limitations  of  the  assumptions  made  in  the  formulation,  offers 
a  powerful  alternative  to  a  purely  experimental  approach  in 
identifying  important  process  mechanisms  and  trends,  which  can 
be  used  to  obtain  better  laser  welds. 


This  paper  examines  the  role  of  traverse  speed,  surface 
reflectivity,  Beer  Lambert  coefficient,  and  changing  liquid 
thermal  conductivity,  on  the  shape  of  the  melt  pool,  and  the 
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cooling  rates  that  occur.  Dependence  of  absorbtion 
coefficient  on  processing  conditions  is  also  examined.  ; 
modified  form  of  a  finite  difference,  heat  conduction 
model[2],  that  now  allows  for  temperature  dependent  substrate 
properties,  in  addition  to  keyholing,  and  surface  convective 
and  radiant  heat  losses  is  used  for  the  purpose.  A  variable 
property  model  is  a  more  realistic  approximation  of  the 
substrate's  response  to  laser  irradiation.  The  model 
calculations  are  of  interest,  as  the  pool  shape  determines  the 
magnitude  and  distribution  of  cooling  rates,  as  well  as  the 
maximum  thickness  that  can  be  welded  in  a  single  pass. 

Cooling  rates,  in  turn,  govern  the  solidification  behaviour, 
and  the  nature,  extent,  and  fineness  of  solid  state  phase 
transformations  that  occur,  including  the  extent  of  the  heat 
affected  zone.  Knowledge  of  the  absorbtion  coefficient  is 
desirable,  for  it  determines  the  laser  power  required  for  a 
desired  melt  depth. 


A  review  of  earlier  models  for  laser  welding  appears  in  [2]. 
A  moving  line-source  model  with  non-uniform  heat  input [31,  and 
a  model  with  heat  flow  in  V  and  U  grooves  in  three-dimensions 
using  variable  properties [4]  are  among  those  proposed.  A 
model  for  surface  melting  with  a  scanning  laser,  that  obtains 
the  temperature  distribution  by  superposing  the  known  Green's 
function  solution  for  the  thermal  distribution  for  a  unit 
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point  source  has  been  reported[5].  Cline  and  Anthony[5] 
assume  constant  substrate  properties,  no  surface  heat  losses, 
and  do  not  account  for  keyholing  or  heats  of  transformation. 
Using  an  approximate  exponential  relation  for  temperature 
penetration  into  the  substrate,  they  estimate  the  depth  of 
penetration  for  laser  welding,  and  report  good  agreement  for 
30U  stainless  steel  for  depths  greater  than  1  cm.  A 
three-dimensional  moving  heat  flux  numerical  model  for  rapid 
melting  and  solidification  of  a  surface  has  recently  been 
publlshed[6].  This  model  provides  information  on  the 
solidification  behaviour  of  laser  melted  materials,  but  Its 
application  to  laser  welding  is  limited  because  keyholing, 
surface  heat  losses  and  variable  substrate  surface  and  bulk 
properties  have  not  been  considered. 


Results  of  present  calculations  show  that  with  an  increase 
in  speed,  the  pool  flattens  out,  is  swept  back,  and  cooling 
rates  increase  near  the  pool  surface  and  edges.  An  increase 
in  Beer  Lambert  coefficient,  fi,  sharply  decreases  the  depth  of 
penetration.  With  the  onset  of  melting,  sensitivity  to 
variations  in  surface  reflectivity  is  reduced.  As  the  thermal 
conductivity  of  the  melt  pool  was  increased,  for  a  given  set 
of  operating  conditions,  shallower  pools,  with  a  higher  aspect 
ratio  were  obtained.  Cooling  rates  increased  as  the  energy 
density  in  the  pool  decreased.  A  dimensionless  melt  front 
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velocity,  0,  i 3  defined,(  see  APPENDIX),  that  provides  an 
upperbound  for  processing  conditions,  such  that  cooling  rates 
exceed  1000  K/s  as  0  approaches  unity.  An  equation  to 
estimate  an  upperbound  for  the  Beer  Lambert  absorbtion 
coefficient,  f},  using  0  is  derived.  The  product  flz,  where  z 
is  the  depth,  was  found  to  be  linear  function  of  the  natural 
logarithm  of  0. 
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3.  THE  MDDIFIED  MODEL 


A  moving  heat  source  model  for  three  dimensional  heat 
conduction  in  the  substrate  during  laser  materials  processing 
is  used.  It  calculates  temperature  fields,  melt  widths  and 
depths  of  penetration.  A  quasi  steady-state  is  assumed  to 
exist,  and  allowance  made  for  keyholing,  laser  beam  mode 
structure,  variable  substrate  properties,  and  surface  heat 
losses  from  a  semi-infinite  slab  of  finite  thickness.  The 
major  modifications  are  the  incorporation  of  variable 
substrate  bulk  and  surface  properties,  and  a  change  in  the 
method  of  solution  to  a  point  successive  overrelaxation 
scheme.  A  linear  grid  is  incorporated  as  an  option,  and  the 
listing  has  been  streamlined  to  reduce  execution  time.  For 
other  details,  see  [2,7]. 


The  assumptions  are: 

1.  The  laser  beam  is  stationary  relative  to  the  earth, 
and  is  incident  at  right  angles  at  the  center  of  the 
substrate  width.  The  substrate  moves  in  the  positive 
x  direction,  (along  the  length),  with  a  uniform 
velocity,  u,  and  is  Infinite  in  length,  but  has 
finite  width  and  thickness. 


35 


8 


I 


2.  Quaai-ateady  3tate  13  assumed. 


3.  The  power  distribution  in  the  beam  is  Gaussian. 


4. 


5. 


The  substrate  thermal  conductivity,  specific  heat, 

thermal  diffusivity  and  the  surface  reflectivity  are 

temperature  dependent.  The  substrate  density  is 

assumed  to  be  independent  of  temperature.  Specific 

heat  variation  accounts  for  heats  of  fusion  and 

vaporization  by  converting  them  into  equivalent  step 

increases  in  values  of  specific  heat,  as  for  example, 

T 

m 


C 


AT 


p  m 


C  dT  +  L 
P 


o 

Surface  reflectivity  is  considered  to  be  zero  when 
the  temperature  exceeds  the  boiling  point. 


1 


6.  There  is  a  radiation  loss  at  the  upper  and  lower 
surfaces  of  the  workpiece. 

7.  Convective  losses  occur  at  the  upper  surface  due  to 
shielding  gas  flow. 

8.  When  the  temperature  at  any  location  exceeds  the 
boiling  point,  it  is  considered  to  have  evaporated, 
but  remains  in  the  conducting  network  at  a 
fictitiously  high  temperature,  to  simulate  the  effect 
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of  high  convection  and  radiation  transfer  effects 
from  the  fast  moving  plasma  in  the  keyhole. 


sdiation  penetrating  the  substrate  is  absorbed 
.^cording  to  the  Beer -Lambert  law, 

P  2  “  Ps  exp(-8z)  2 

where  $  is  the  absorbtion  coefficient  (m  )  and  is 
considered  independent  of  position  within  the 
keyhole,  and  Pz  is  the  absorbed  power  density  at 
depth  z  which  was  of  value  Ps  at  the  surface. 
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4.  MODEL  CALCULATIONS 


Calculations  were  made  for  a  Ti-6AL-4V  substrate.  Data  on 
property  variation  (Fig.  1-3)  in  the  solid  3tate  are  from 
published  1 iterature[S ] ;  projections  above  melting  point  are 
thought  to  be  reasonable  approximations  of  the  expected 
behaviour.  Actual  data  for  these  high  temperatures  are 
unavailable. 


Changes  in  pool  shape  and  variations  in  cooling  rates  were 
the  two  areas  of  interest.  Changes  in  depth  of  penetration, 
and  in  surface  melt  widths  were  used  to  quantify  variations  in 
pool  shape.  Computer  calculations  were  made  to  find  the 
effects  of  changes  in  traverse  speed.  Beer  Lambert 
coefficient,  effective  thermal  conductivity  in  the  liquid 
state,  and  surface  reflectivity  upon  pool  shape.  Overall 
rates  of  heat  transfer  are  increased  by  convective  heat 
transport,  due  to  fluid  flow  in  the  melt  pool.  This  effect  is 
approximated  by  arbitrarily  increasing  the  effective  thermal 
conductivity  in  the  liquid  state  in  the  model  as  it  assumes  a 
static  pool.  Temperature  gradients  in  the  x-direction,  and 
the  traverse  speed,  were  used  to  determine  the  variations  in 
cooling  rate  with  melt  depth  and  traverse  speed.  The  cooling 
rates  were  related  to  the  energy  density  in  the  pool,  assuming 
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good  coupling  between  the  laser  and  substrate.  A 
dimensionless  melt  front  velocity,  ^>,  was  defined,  where 

$=  (P/d2)/(upC  AT  )  3 

*■  pm 

Using  it  is  possible  to  compute  an  upper  limit  to  the  value 
of  the  absorption  coefficient,  (see  APPENDIX),  and  effective 
limiting  values  for  different  energy  densities  in  the  pool 
are  calculated.  The  value  of  the  depth  of  penetration  used  to 
calculate  $  here,  are  derived  from  the  model  itself,  for  a 
constant  £  value  of  1000. 


i 

It  should  be  noted,  that  P/ud  is  only  a  measure  of  the 
actual  energy  density  in  the  melt  pool,  as  it  assumes  that  the 
melt  pool  is  a  rectangular  parallelopiped,  of  crossectlon  d1, 
of  length  u,  which  absorbs  W  joules  of  energy  every  second. 

It  assumes  that  no  losses  occur  into  the  substrate,  and  that 
laser'-substrate  coupling  is  perfect.  In  practice,  under 
steady  state  conditions,  the  actual  energy  density  will  be  a 
constant  fraction  of  this  estimate  (see  also  [93). 
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5.  MODEL  RESULTS 


With  an  increase  In  speed,  the  the  pool  flattens  out,  is 
swept  back,  and  the  leading  edge  makes  a  smaller  angle  with 
the  center  line  (Fig.  4-6}.  The  point  of  highest  temperature 
shifts  behind  the  laser  beam.  The  width  decreases  too,  but 
not  as  sharply  as  the  depth.  For  a  given  incident  power 
density,  an  increase  in  speed,  increases  cooling  rates,  with 
higher  rates  being  observed  near  the  surface  and  the  edges  of 
the  melt  pool  (Fig.  7,8).  An  increase  in  f}  results  in  sharp 
decrease  in  depth  (Fig.  9),  while  the  width  of  the  pool  is 
not  affected  as  much.  If  the  surface  reflectivity  is  too 
high,  (0.99  say),  no  melting  occurs  (Fig.  10).  However,  if 
the  reflectivity  is  such  that  the  fraction  of  energy  absorbed 
is  sufficient  to  initiate  melting,  then,  the  sensitivity  to 
changes  in  reflectivity  is  greatly  reduced. 


Increased  effective  thermal  conductivity  in  the  liquid  state 
can  significantly  affect  the  shape  of  the  weld  puddle  (Fig. 
11).  The  aspect  ratio  of  the  pool  increases  as  the  Peclet 
number  (ud/ct  )  decreases.  The  decrease  in  depth  implies  that 
materials  with  a  higher  thermal  conductivity  are  harder  to 
weld,  as  is  observed  in  practice,  for  a  given  set  of  operating 
conditions . 
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As  the  energy  density  (J/nm*  )  of  the  pool  decreases,  cooling 
rates  increase  (Fig.  12).  Cooling  rates  in  excess  of  ten 
thousand  degrees  a  second  are  obtained  when  values  of  the 
dimensionles  melt  front  velocity  are  less  than  a  hundred  (Fig. 
13) •  Changes  in  energy  density  in  the  pool  also  affect  the 
values  of  the  Beer  Lambert  absorbtion  coefficient,  /S,  with  the 
product  of  fi  and  z,  the  depth  of  penetration,  being  linearly 
related  to  the  natural  logritho  of  <P  ,  the  dimensionless 
energy  density,  (Fig.  14). 
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6.  DISCUSSION 


The  paucity  of  high  temperature  material  property  data  is  a 
source  of  major  difficulty  in  using  variable  property  models. 
An  assumption  of  constant  properties  is  unrealistic,  and  leads 
to  predictions  of  doubtful  value.  In  the  present  model,  the 
variation  in  thermal  conductivity  at  high  temperature  reflects 
the  expected  effects  of  convective  heat  transport  due  to  fluid 
flow  in  the  melt  pool,  and  the  rapid  heat  transfer  in  the 
plasma  filled  keyhole.  Specific  heat  variation  accounts  for 
heats  of  fusion  and  vaporization.  The  heats  of  transformation 
are  easily  converted  to  specific  heat  equivalents.  The 
thermal  diffusivity  is  derived  from  the  conductivity  and 
specific  heat,  assuming  a  constant  substrate  density.  While 
the  arbitrary  choice  of  these  temperature  dependencies  is 
unavoidable,  the  assumed  temperature  dependence  is  crucial  to 
model  results.  Hence,  it  is  the  trend  that  emerges,  rather 
than  the  exact  numerical  value,  which  is  of  greater  importance 
here. 


As  the  traverse  speed  lncre"  ',  the  time  of  interaction 
between  the  laser  and  the  substrate  decreases.  This  results 
in  shallower,  colder  pools,  as  the  substrate  absorbs  a  smaller 
amount  of  the  Incident  energy.  Also,  front-to-back  velocity 
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of  the  molten  material  increases,  so  that  convective  heat 
transport  past  the  beam  increases,  resulting  in  pools  that  are 
swept  bade.  Both  these  factors  lead  to  Fig.  4-6.  The 
shallower,  colder  pools  have  a  higher  surface  area  per  unit 
volume,  in  contact  with  a  relatively  cooler  substrate.  They 
are  thus  subjected  to  higher  rates  of  heat  extraction  by  the 
substrate,  and  cooling  rates  increase.  In  general ,  cooling 
rates  Increase  as  solidification  proceeds  towards  the 
surface[10],  as  happens  at  the  trailing  edge  of  the  pool. 

While  this  Information  is  not  entirely  new(see  [5], [11]),  it 
does  serve  as  a  check  on  model  predictions. 


Laser  radiation  couples  to  a  substrate  better  when  it  is 
molten  rather  than  solid[12].  Thus,  small  variations  in  solid 
reflectivity  will  have  a  lesser  effect  if  the  laser  is 
incident  upon  a  pool  of  molten  material  at  the  surface  (Fig. 
10).  &irface  reflectivity  variations  will  be  more  critical 
during  laser  heat  treatment,  when  no  surface  melting  occurs, 
than  during  laser  welding  or  alloying. 


Since  this  model  does  not  explicitly  allow  for  fluid  flow, 
the  increased  convective  heat  transport  away  from  the  beam  can 
only  be  acoounted  for  in  an  approximate  way,  by  increasing  the 
effective  liquid  thermal  conductivity(Fig.  11).  The  trend 
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there  13  reasonable ,  In  that  the  aspect  ratio  increases  as 
fluid  flow  becomes  more  vigorous,  or  effective  thermal 
conductivity  increases.  Actually,  flow  at  the  surface  will 
increase  melt  widths,  and  reduce  depths,  as  more  material  is 
melted  by  heat  transported  laterally  away  from  the  beam.  It 
also  implies,  that  if  an  arbitrary  increase  in  effective 
liquid  conductivity  is  assumed  in  models  with  constant 
substrate  properties [6],  then  the  choice  largely  dictates  the 
predicted  melt  pool  shapes. 

A  measure  of  the  rate  at  which  the  solid-liquid  (S-L) 
interface  sweeps  through  the  material  is  provided  by  the 
cooling  rate.  In  laser  welding,  solidification  occurs  at  the 
S-L  interface  trailing  the  beam,  and  moving  with  it.  The 
cooling  rate,  except  for  starting  and  stopping  transients,  is 
determined  by  the  rate  at  which  this  part  of  the  S-L  interface 
travels  through  the  substrate.  As  <J>  approaches  1,  the  S-L 
interface  will  sweep  through  the  material  at  faster 
velocities,  increasing  cooling  rates.  Hence,  ^  may  be  used  to 
estimate  a  working  range  of  processing  conditions,  where  we 
can  expect  high  cooling  rates  (Fig.  13). 


0  is  also  the  ratio  of  the  energy  density  in  the  pool  to 
that  required  to  cause  melting.  As  *  approaches  1,  the 
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superheat  in  the  pool,  decreases,  so  that  the  pool  will  be 
subjected  to  higher  rates  of  heat  extraction,  Increasing 
cooling  rates(Fig.  12). 


In  modeling  of  laser  processes,  as  well  as  in  practice,  a 
knowledge  of  the  absorbtion  coefficient,  /?,  is  desirable.  In 
as  much  as  has  not  been  directly  measured,  the  equation  for 
estimating  an  upperbound  on  £  should  be  quite  useful.  A 
better  estimate  can  be  obtained,  if  actual  depths  of 
penetration  are  used,  with  the  energy  density  required  to 
initiate  melting.  Both  these  quantities  are  readily  measured 
experimentally.  $  can  be  estimated  if  the  depth  of 
penetration  is  known,  and  the  constant  A  is  estimated  from 
experiments,  by  using  the  relationship  plotted  in  Fig.  14. 


The  results  of  this  model  are  relevant  to  laser  welding. 

The  effects  of  a  change  in  the  absorbtion  coefficient,  its 
variation  with  pool  energy  density  provide  new  insights.  A 
means  to  estimate  an  upperbound  for  (&  is  useful.  The  role  of 
varying  reflectivity,  and  increased  fluid  thermal  conductivity 
in  the  pool  point  to  important  mechanisms.  A  knowledge  of  the 
variation  of  cooling  rates  with  energy  density,  and  a  means  of 
selecting  operating  conditions  that  will  produce  high  cooling 
rates  via  <P  should  help  in  the  production  of  better  welds. 
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7.  CONCLUSIONS 


The  following  conclusions  can  be  drawn  from  the  above 
results: 

1.  The  absorbtion  of  laser  energy  in  the  keyhole,  as 
determined  by  the  Beer  Lambert  coefficient  p,  affects 
the  depth  of  penetration  more  than  the  width  of  the 
weld. 

2.  Variations  in  the  surface  reflectivity  are  less 
critical  in  laser  welding  or  alloying,  as  surface 
melting  occurs. 

3.  The  aspect  ratio  of  the  weld  pool  increases,  as 
effective  liquid  thermal  conductivity  increases,  with 
an  increase  in  the  vigor  of  fluid  flow  in  the  melt 
pool. 

4.  Cooling  rates  Increase  as  the  energy  density  in  the 
melt  pool  decreases. 

5.  An  operating  region  of  high  cooling  rates  may  be 
estimated  using  the  dimensionless  melt  front 
velocity,  0000,  where 

(p/d^)/(upc  AT  )  4 

P  a' 
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6.  An  upperbound  for  d  may  be  estimated,  using 

a  =  A  (A  i)  5 

z 

where  the  constant  A  can  be  determined  from  the 
experimentally  determined  power  density  required  to 
initiate  melting  at  that  speed  as 
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portion  of  the  substrate  itself,  melting  will  occur  only  if 

the  irradiation  intensity  is  greater  than  the  energy  density 

required  to  produce  melting.  The  energy  density  at  the  bottom 

of  the  pool  is  AP  /u,  where  A  is  a  constant,  u  13  the  traverse 
z 

2 

speed,  m/s,  and  ^  is  the  power  density,  W/m  .  For  melting 

to  occur,  AP  /u  must  at  least  equal  pH  AT  .  ■  Using  the 

z  Pm 

Beer  Lambert  law,  and  substituing  for  P  ,  we  obtain 


u  p  CpAT^  ^  exp(-Bz) 


A'4 


Recognizing  that  Ps/upCpATm  is  just  the  dimensionless  melt 
front  velocity,  4,  we  obtain 


Bz  «*  inCA© 

The  constant  A  may  be  estimated  if  the  power  density  to 

2 

initiate  melting, (p/d  )  at  that  speed  is  determined 

expt 

■experimentally.  Then, 


A'S 


A  -  (upCpATm)/(P/d2)expt  A-6 

In  the  keyhole,  the  reflectivity  is  assumed  to  be  zero,  due  to 
internal  reflections.  In  actual  practice,  there  will  be  a 
layer  of  liquid  between  the  bottom  of  the  keyhole  and  the 
unmelted  substrate.  Thus,  excess  heat  energy,  will  have  to  be 
supplied  at  the  bottom  of  the  keyhole,  to  heat  fresh  substrate 
through  the  liquid  layer  at  the  floor  of  the  keyhole.  Hence, 
this  estimate  of  B  is  an  upperbound. 
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10.  LIST  OF  CAPTIONS 

1.  Variation  with  temperature  of  thermal  conductivity 
for  T1-6A1-4V.  Data  above  melting  point  are 
arbitrary. 

i 

2.  Variation  with  temperature  of  specific  heat  for 
T1-6A1-4V.  Data  above  melting  are  arbitrary,  and 
have  been  adjusted  for  heats  of  melting  and 
vaporization. 

3.  Variation  with  temperature  of  thermal  diffusivlty  for 
T1-6A1-4V,  derived  from  Fig.  1  and  2,  assuming 
constant  substrate  density. 

4.  Predicted  melt  pools,  assuming  Beer  Lambert 
coefficient  of  1000  m'1. 

5.  Predicted  melt  pools,  assuming  Beer  Lambert 
coefficient  of  1000  m'1  . 

6.  Predicted  melt  pools,  assuming  Beer  Lambert 
coefficient  of  1000  m'1  . 

7.  Center  line  cooling  rate  at  the  surface,  for  laser 
power  of  2k W,  and  beam  diameter  of  0.2  mm.  Note  how 
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point  of  peak  temperature  shifts  with  increasing 
apeed.  Cooling  rates  plotted  for  trailing  half  of 
pool.  Abcissa  3hows  distance  behind  beam  center. 

8.  Center  line  cooling  rate  0.4  mm  below  surface,  for 
laser  power  of  2kW,  and  beam  diameter  of  0.2  mm. 

Note  how  point  of  peak  temperature  shifts  with 
Increasing  speed.  Cooling  rates  plotted  for  trailing 
half  of  pool.  Abcissa  shows  distance  behind  beam 
center . 

9.  Predicted  melt  fronts  for  different  values  of  Beer 
Lambert  coeffient. 

10.  Predicted  melt  fronts  for  different  values  of  surface 
reflectivity.  Note  that  for  reflectivity  of  0.99, 
model  predicts  no  melting  will  occur.  Assumed  of 
Beer  Lambert  coefficient  1000  m"' . 

11.  Predicted  melt  fronts  for  different  values  of 
effective  thermal  conductivity. 

12.  Variation  of  center  line  cooling  rate  near  beam,  with 
estimated  melt  pool  energy  density. 

13«  Variation  of  center  line  cooling  rate  near  beam,  with 


52 


25 


dimens tonless  melt  front  velocity. 


14.  Variation  in  product  of  Beer  Lambert  coefficient,  , 
m'1  ,  and  z  the  depth  of  penetration,  m,  with  changes 
in  dimensionless  energy  density, ,  for  different 


values  of  A. 
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11.  LIST  OF  SYMBOLS 


A  Constant  relating  measured  and  calculated 

energy  density  of  melting 

Cp  Specific  heat  of  substrate 

C  Specific  heat  of  substrate  corrected  for 
P 

heat  of  fusion 
d  Beam  diameter 

k  Thermal  conductivity  of  substrate 

L  Latent  heat  of  fusion  of  substrate 

P  Incident  laser  power 

Ps  Incident  energy  flux  at  surface 

Pz  Absorbed  power  density  at  depth  z 

T  Temperature 

Tq  Ambient  temperature 

T^j  Melting  point 

AT  (T  -T  ) 

m  o 

u  Traverse  speed 

z  Melt  depth 

a  Thermal  diffusivity  of  substrate 
0  Beer  Lambert  absorbtion  coefficient 
P  Substrate  density 

^  Melt  front  velocity  or  pool  energy  density 


no  unit 


J/kg  K 
J/kg  K 


U/m  K 

J/kg 

W 

W/m2 

W/m2 

K 

K 

K 

K 

m/s 

m 

mVs 

.i 

m 

kg/m3 
no  unit 
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Variation  with  temperature  of  thermal  conductivity 
for  T1-6A1-IIV.  Data  above  melting  point  are 
arbitrary. 


Temperature  Fit 


Variation  with  temperature  of  thermal  dlffusivlty  for 
Tl-6«-l|V,  derived  from  FI g.  1  and  2,  assuming 
constant  substrate  density. 


Predicted  Belt  poola,  aaaualng  Beer  Leabert 
coefficient  of  1000  a"'. 


_ £  4 


Predicted  Belt  pools,  essuslng  Beer  Leabert 


2.0x10 
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Center  line  cooling  rate  at  the  surface,  for  laser 
Vow"  of  ZcV ,  and  beaa  diameter  or  0.2  mo.  Mote  how 
Point  of  peak  temperature  shifts  with  Increasing 
»pee<3.  Cooling  rates  plotted  for  trailing  half  of 
P°ol.  Abcissa  shows  distance  behind  beam  center. 
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1*1.  Variation  In  product  of  Beer  Laabert  coefficient,^ 
a  ,  and  x  the  depth  or  penetration,  m,  with  changes 
In  dlnenslonless  energy  density,  for  different 


A  .  .'»rm rp,  *' . . i  toer. ,  P.  henry ,  G.C.  Lim 
Imperial  'ol  lege, London , England . 


INTRODUCTION 


H i gh  power  Inner  welding  is  now  an  accepted  material  processing 
■  nn  i '}ue.  In  gaining  this  acceptance  a  large  amount  of  research  has 

been  done  on  the  effect  of  various  parameters  e.  g.  traverse 

spool, laser  power, etc.  /I, 2/. One  area  which  has  received  little 

attention  is  the  effect  >f  mode  structure  if  the  laser  beam  on  the 

wold  profile. 

There  are  two  main  reasons  why  sucn  an  investigation  is  import¬ 
ant.  First,  many  of  the  larger  power  machines  do  not  produce  a 

Gaussian  beam, and  so  the  significance  of  a  non-Gaussian  beam  is 
important  in  assessing  the  use  of  these  machines.  Secondly , since  the 
mode  structure  of  the  output  from  a  machine  may  vary  over  a  period  of 
time.it  is  necessary  to  have  a  knowledge  of  the  effects  of  the  mode 
structure  so  that  the  reproducibility  of  the  weld  can  be  guaranteed. 

As  a  starting  point  of  this  investigation  , the  profile  of  a 
focussed  beam  operating  in  a  TEMmn  mode  is  shown  to  have  the  same 
profile  as  the  unfocussed  beam  when  aperture  effects  are  neglected  and 
perfect  optics  are  assumed.  From  this  it  can  be  deduced  tnat  the 
focussed  beam  has  the  same  intensity  as  the  raw  beam  for  any  intensity 
profile. 

A  new  definition  of  the  spot  size  is  introduced  which  allows  a 

consistent  comparison  of  the  spot  sizes  of  the  different  modes.  Using 

a  finite  difference  numerical  model  of  the  laser  welding  process  an 
investigation  is  made  into  whether  effect  of  the  mode  structure  on 
the  weld  profile  is  due  soley  to  the  varying  effective  spot  sizes  of 
the  modes, or  whether  the  actual  mode  structure  itself  has  an  effect. 

A  theoretical  investigation  is  made  into  the  effect  of  a  finite 
lens  aperture  on  the  spot  size, and  it  is  also  shown  how  lens 
aberrations  can  be  taken  into  account.  Finally  an  experiment  is 
described  to  validate  the  theoretical  results  outlined  in  this  paper. 

FOCUSSING  OF  LASER  BEAMS 


The  TEMmn  modes  for  a  laser  with  cylindrical  symmetry  are  given 
in  /V  as 


O) 


where  (r,9)  are  polar  coordina tea , w  is  the  characteristic  radius  of 
the  beam.  L  is  the  Laguerre  polynomial  of  order  n  and  degree  m,  and 
Ann  is  the  amplitude  of  the  TEMmn  mode. The  intensity  is  given  by 


K  e.e )  =  I/O1 


in  1  i  :*  1 
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where  th»?  geometry  jrivolv*#a  il 
focal  length  of  th«-  Ions,. ml  ■ i  is 


... 

,  it 


i(i‘ -  ,>ir'(l"',))  , 

6  *  f  arc 

i j - *  fit*-  1  ui  rig.l.  k  =  ,’ir/X 
t,n<-  aperture  ri'liuj. 


,  f  is  the 


Qm»(M)  =  e  L^s  ^X(zr)rAf 


(4) 


f:C 

whore  z=kp/s.  At  foe  as  the  complex  exponential  term  disappears. 
Aperture  effects  can  be  neglected  by  taking  a  as  infinite.  In  this 
case  it  can  be  shown  that 


(3mn  (f,^)  =  Amn  (  ^  ,$) 


where 


w. 


=  M 


»rw 


is  the  characteristic  radius  of  the  focussed  beam.  Since  the  Amn  form 
an  orthogonal  set  of  functions  and  (3)  defines  a  linear  integral 
operator  on  the  functions  Amn, then  it  can  be  shown  that 


a(e>*)  - 


(6) 


for  ali  beams  with  unfocussed  amplitude  A(r,6) 


SPOT  SIZE 


There  have  been  many  definitions  introduced  to  describe  the  spot 
size  of  the  laser  beam.  The  characteristic  radius, w,is  not  a  good 
measure  of  the  spot  size  since  the  effective  spot  size  of  the  higher 
order  modes  is  much  greater  than  the  characteristic  radius.  The 
important  point  about  the  characteristic  radius  is  that  it  is  a 
function  of  the  geometry  of  the  laser  cavity  alone  /3/,  and  is  there¬ 
fore  independent  of  the  mode  structure. 

The  one  drawback  with  the  characteristic  radius  and  other 
definitions  is  that  the  power  contained  within  the  defined  radius  is 
dependent  on  the  mode  itself. 

To  rectify  this  the  following  definition  of  spot  size  is  intro¬ 
duced.  A  laser  beam  of  total  power  P  and  intensity  distribution 
l(r,C)  has  radius, R,  defined  by 

J*  f  Ur,e)rdrAb  =  U-e'^P  (7) 

r  *c  fc-v 

’•i"  ?  i .  -  *  >  r  (!-•“*'  ;«»  f  •>.■»♦  t  *i«*  riiiu.'.  » s  »*  pi  t 
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Mod"  factors  for  the-  lower  order  modes  are  given  in  table  1  below 


TABLE  1 

Mode  factors  of  TEMmn  modes 


—  N  — 
0  1 


1 

0 

1  .00 

1  .32 

1  -56 

I 

M 

i 

1 

1  .6b 

1  .86 

2.08 

1 

2 

2.12 

2. Ml 

2 . 48 

APERTURE  EFFECTS  AND 

ABERRATIONS 

It  has  been  shown  that  the  focussed  beam  profile  is  that  of  the 
raw  beam  if  aperture  effects  are  neglected , and  so  the  spot  size  will 
be  given  by  the  characteristic  radius  multiplied  by  the  mode  factor. 
The  integral  in  (4)  was  calculated  numerically  to  show  the  effect  of  a 
finite  aperture  on  the  spot  size  for  the  low  order  modes.  The  results 
are  shown  in  table  2. 


TABLE  2 

Increase  in  spot  size  due  to  finite  aperture 


a/w 

m=0, n=0 

m=0, n=1 

m=1 ,n=0 

m=  1  , n= 1 

0.5 

2.55 

2.87 

1  .64 

2.04 

1  .0 

1.37 

1 .61 

1.55 

1.28 

2.0 

1  .00 

1.00 

1.00 

1.00 

It  can  be  seen  that  if  the  aperture  radius  is  greater  than  twice 
the  characteristic  radius  of  the  beam  then  the  effect  is  negligible. 


Aberrations  can  be  accounted  by  including  an  extra  factor 

exp(ikW(r,e))  (9) 

in  (4),W  is  the  aberration  function  /5/  and  describes  the  deviation  of 
the  actual  wavefront  from  a  spherical  wavefront.  For  example! 

W(r,6)=Cr*  (io) 

describes  spherical  aberration.  C  depends  on  the  amount  of  aberration 
present.  Since  the  beam  extends  over  a  wider  range  for  higher  order 
modes.it  would  be  expected  that  higher  order  modes  would  be  affected 
more  by  aberrations. 

THE  EFFECTS  OF  MODE  STRUCTURE  ON  LASER  WELDING 

Having  detailed  the  relationship  between  the  raw  beam  ind  the 
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In  order  to  see  if 

the 

actual  inode 

structure  itself  has  an  effect 

the  model  was  run  for  a 

TEMOt*  mode  and 

for  a 

Gaussian  mode  of  the 

same  effective  radius 

,as 

well  as  a 

Gaussian  beam  of  the  same 

characteristic  radius, table 

3 

TABLE  3 

Effects  of  mode  structure 

jn  weld 

profile 

(a 

1 1  lengths  in 

mm ) 

P  = 

1600W 

2000V 

w= 

0.25 

0.  50 

w 

d 

w 

d 

TEM01*  w 

o 

CO 

1  .7 

2.0 

1.4 

TU100  w 

- 

- 

1  .3 

2.0 

Traoo  1.32W 

o 

CO 

2.2 

2.3 

2.0 

It  can  be  seen  that  the  TEM01*  beam  has  less  penetration  than  a 
Gaussian  of  the  same  effective  radius.  This  suggests  that  the  mode 
structure  does  have  an  effect  apart  from  just  its  increase  in  size, but 
it  is  felt  that  further  investigation  is  necessary. 

EXPERIMENT 


Experiments  are  in  progress  at  Imperial  College  to  validate  the 
ideas  expressed  in  this  paper. The  experiment  simply  consists  of 
placing  an  ALL  Laser  Beam  Analyser  in  the  path  of  a  beam  from  a 
Control  Laser  Ltd.  2kW  CO  laser, so  that  the  mode  structure  of  the  raw 
beam  can  be  monitored  while  welding  is  carried  out. 

The  ALL  beam  analyser  can  also  be  used  to  measure  the  focussed 
beam  profile, and  this  will  be  used  to  check  the  relationship  between 
the  raw  and  focussed  beams. 
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I  A;  I)  V-  NSIONAl  TRANSIENT  MODEL  FOR  CONVECTION 
'f.  ASP  M  VP  LTFD  POOLS 

•.,  *  ,  V|/  •  l****,  jnd  W  M  Chfo’" 


A  t*3 -aine/is i ona  i  transient  model  fo r  convective 
neat  transfer  and  surface  tension  driven  fluid  flew  Is 
developed.  The  model  describes  the  transient  behavior 
of  the  heat  transfer  process  of  a  stationary  band 
source.  Some  semi -quant i tat i ve  understanding  of 
scanning  is  obtained  by  a  co-ordinate  transformatl on. 
The  non-dimens  Iona  1  form  of  the  equations  are  derived 
and  four  dimensionless  parameters  are  identified, 
namely,  Pec  let  number  ( Pe ) ,  Prandtl  number  (Pr) , 
Surface  tension  number  (S),  and  dimensionless  melting 
temperature  (T*  .  ).  Their  governing  characteristics 
and  their  eff^it  on  pool  shape,  cooling  rate,  velocity 
field  and  solute  redi stri buti on  is  discussed.  Numerical 
solution  is  obtained  and  presented. 

NOMENCLATURE 

C  constant  which  defines  the  interface, 

no  unit,  dimensionless 

d  diameter  of  laser  beam;  mm,  [L] 

D  ff  effective  diffusion  coefficient; 

mz/sec.  [LZ/t] 

k  thermal  conductivity;  kW/m*K,  [Ml/tZT] 


rQ  radius  of  laser  beam;  mm,  [l] 

S  Surface  tension  number;  no  unit,  dimensi onless 

T  temperature;  °K,  [T] 

Tmelt  melting  temperature;  °K,  [T] 

T*  Olmensionless  melting  temperature 

melt  r  -  T 

(  metd1  );  no  unit,  dimensionless 

qd  /  K 

T  temperature  of  metal  when  It  Is  not  heated; 

met4’  x  m 

u  vetodty  vector;  mm/sec,  [l/t] 

u  ii-component  of  u;  mm/sec,  [l/t] 

u  scanning  speed  of  the  laser  beam; 

mm/sec,  [l/t] 

v  y-component  of  u;  mm/sec,  [l/t] 

w  z-component  of  u;  mm/sec,  [l/t] 

«,y,z  Cartesian  coordinate;  mm,  [l] 


t  length  of  the  rectangular  heat  source;  mm,  [l]  Grdek 

p  pressure;  N/m2.  [M/t2l]  »  density;  kg/m3.  [M/l3] 

Pe  Pec  let  number  u  d/<,  no  unit,  dimensionless  *  thermal  diffusivity,  m  /sec.  [L  /t] 

o 

q  net  heat  flu«  from  laser;  kw/mm2,  [M/t2]  «  kinematic  viscosity;  m2/sec.  [l2/t] 

u  viscosity;  kg/sec-*.  [M/t-l] 

cj  surface  tension,  N/m,  [M/tZ] 


Re  Reynolds  number  uQ  d/v;  no  unit,  dimensionless 
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pool  and  the  melt  pool  encountered  in  laser  surface 
alloying. 

In  this  paper,  the  convective  heat  transfer  and 
fluid  flow  in  a  laser  melted  pool  is  analyzed.  ft 
transient  two-dimensional  model  for  surface  tension 
driven  fluid  flow  in  laser  melted  pool  was  developed. 
The  rode!  describes  the  transient  behavior  of  the  heat 
transfer  process  of  a  stationary  band  source.  Some 
semi -quantitative  understanding  of  scanning  can  be 
obtained  by  a  coordinate  transformation.  Four 
dimensionless  parameters  arising  from  the  general 
three-dimensional  mathematical  formulation  are  scanning 
Peclet  number  (Pe),  Prandtl  number  (Pr),  Surface 
tension  number  (S).  and  the  dimensionless  melting 
temperature  (T*  .  ).  Each  of  these  would  govern  the 
characteristic's8^  the  problem  and  each  has  its  <xn 
physical  interpretation,  numerical  solutions  using  the 
program  SOLA  [9]  are  obtained  for  five  cases  and  their 
results  are  presented.  The  recirculation  fluid  flow 
can  be  observed  clearly.  The  prediction  of  the  pool 
shape,  in  terms  of  aspect  ratio  (width  divided  by 
depth),  as  a  function  of  the  governing  parameters,  Pr 
and  S,  are  plotted  and  discussed.  The  cooling  rate 
along  the  interface  is  calculated  and  its  implication 
to  microstructure  is  discussed.  Some  semi -quantitative 
prediction  of  mixing  is  also  discussed. 

PHTSiCS  OF  THE  PROCESS 

ft  two-dimensional  analysis  is  presented,  modeling 
a  band  source  laser  beam  having  a  defined  power 
distribution  striking  the  surface  of  an  opaque  material 
of  infinite  width,  thickness,  and  length  (Fig.  1). 
Some  of  the  incident  radiation  is  reflected  while  the 
rest  is  absorbed.  The  heat  absorbed  developed  a  molten 
pool.  Owing  to  the  high  temperature  gradient  on  the 
surface  in  the  transverse  direction,  the  surface 
tension  gradient  is  developed.  It  is  this  mechanism 
that  drives  the  flow.  As  the  flow  develops,  energy 
transfer  nechanism  becomes  a  convection  problem  with 
the  flow  driven  by  the  surface  tension  gradient.  The 
basic  assumptions  are: 


1.  ’m*  . :  j  t  *v  t  i  '•  *  s  cons  1  dered  to  >-■> 

.ha?  t h,-.  Navi  or- StoVes  equation  is  ap,::  i  ahtp. 

2.  All  the  properties  of  the  liquid  metal  jn.l  s  >  f ' 

metal  are  constant,  independent  o!  temperature 
(except  the  surface  tension).  This  allots 

simplifications  of  the  model;  however,  variable 
properties  can  he  treated  with  slight  modifica¬ 
tions.  The  dependence  of  surface  tension  on 

tempera  tore ,  the  driving  force  of  the  flow,  is 
assumed  to  be  linear. 

3.  Laton  heat  of  fusion  is  neglected  since  tne  energy 

liberated  is  small  compared  to  total  enthalpy 

change  associated  with  temperature  differences. 

4.  Thermal  conductivity  is  assumed  to  be  the  same  for 
both  liquid  and  solid  phases  for  the  save  of 
simplicity  of  the  model. 

5.  Surface  of  the  melt  pool  is  assumed  to  be* flat,  to 

simplify  the  surface  boundary  condition  and, 

hence,  surface  rippling  is  neglected. 

MATHEMATICAL  FORMULATION 

The  appropriate  governing  equations  are: 

Energy  equation, 

ST/Jt.  .  (u  •  V)  T  =  <  V*  T  (I  ) 

Continuity, 

V  -  u  =  9  (2) 

Momentum  equations, 

3u/H  ♦  (u  •  V)  u  *  -( l /p) ( 3P/3x)  ♦  v  u  (3) 

aV/dt  ♦  (u  •  7)  V  *  -(l/p)(3P/3y)  ♦  v  V*  V  (4) 

and  boundary  conditions  are  at  the  surface 
y  =  0,  v  =  0,  u(IW3y)  *  -o’(3T/3x) 

1  ro!  otberv)  se 

and  at  tbe  liquid-solid  interface 
f (»,y)  -  constant;  u  *  v  «  U.  T  *  (6) 
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Peclet  number  can  be  interpreted  as  the  ratio  of 
heat  diffusion  time  to  scanning  interaction  time. 
Reynolds  number  is  the  quotient  of  Reynolds  number  and 
Prandtl  number.  The  latter  is  a  property  of  the  molten 
material --the  ratio  of  momentum  diffusion  and  the 


me 

problem  becomes  conduction.  On  the  other  hand  if  the 
number  is  very  large,  convection  flow  would  be  very 
vigorous  and  the  physics  of  the  problem  is  quite 
different.  In  fact,  it  becomes  convection  dominated. 
The  shape  of  the  molten  region  tends  to  be  wider. 
Dimensionless  siting  temperature  would  govern  the  sue 
of  the  molten  region.  The  higher  the  number,  smaller 
the  size. 


Table  1  Physical  Property  of  Steel,  Al ,  and  NaNC^ 


Steel 

Al 

NaNo3 

p 

7.0  x  103  kg/m3 

2.38b  x  103  kg/m3 

1.904  x  103  kg/m3 

u 

5.6  *  10'3  N  s  m'2 

2.8  x  10"3  N  s  m_? 

*2.78  x  I0"3  N  s  n" 

V 

7.84  x  I0'6  m2/sec 

1.17  x  10*^  m^/sec 

1.46  x  10"®  m2/sec 

k 

31  W/n°i: 

100  W/m°  K 

1  .  0  x  10*®  W/m’K 

K 

1 .0  x  10"®  m2/sec 

6.0  x  10“^  m^/sec 

5.6  x  10"7  m2/sec 

3melt 

1800°C 

600°C 

306. 8°C 

0* 

-0.112  x  10'3  N/m/°K 

-0.3S  x  10"3  N/m/°K 

-0.7  x  I0'4  N/m/°K 

Ta t>  1  e  2  Order  of 

Magnitude  of  Parameters 

Steel 

Al 

Ha  Mo  j 

Rt 

2.55  to  8. 29 

17  to  55 

14  to  45 

Re 

2  to  6.5 

0.3  to  1.0 

36  to  120 

S 

5  x  I04  to  2.5  i  I05 

10S  to  5  x  10® 

2.0  x  I08  to  I09 

TAeH 

0.005  to  0.001 

0.006  to  0.016 

3  x  10"10  to  7  x  !0"; 
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a.  Surf  j c e  Tempera*.  u«*e 

T ht?  Suffice  teiverdt jre  within  t nu  molten  region, 
for  different  surface  tension  number*  and  Prandtl 
Oj-oers  are  plotted  m  Fig.  8  and  9  respectively.  The 
curve  e  x  n  i  b  1 1  $  the  usual  feature ,  >  t  attends  1 1  s 
mauling™  at  the  center  and  decreases  as  it  moves  to  the 
edge.  The  temperature  gradient,  driving  force  of  the 
fluid  flow,  is  greatest  at  the  edge  of  the  laser  beam 
{x*  =  O.S).  The  effect  of  surface  tension  number  on 
the  surface  temperature,  as  can  be  seen  in  Fig.  8,  is 
rather  small  as  compared  to  the  effect  of  Prandtl 
number  on  it.  Fig.  9  Prandtl  number,  as  explained 
earlier,  is  the  ratio  of  momentum  and  thermal 
diffusivity.  An  increase  in  Prandtl  number  can  be 
interpretted  as  a  decrease  in  thermal  diffusivity. 
Thus,  the  material  is  heated  up  more  slowly  when  its 
Prandtl  number  is  higher.  As  can  be  seen  in  Fig.  9, 
the  temperature  gradient  decreases  with  Prandtl 
number.  This  causes  the  magnitude  of  the  recirculating 
flow  decrease  with  Prandtl  number. 

C.  Horizontal  Surface  Velocity 

The  horizontal  surface  velocity  at  t*  »  1  and 
x*  =  O.S  as  a  function  of  S  number  is  plotted  in 
Fig.  10.  The  nvjgnitude  increases  with  S  number  which 
is  expected.  The  horizontal  surface  velocity  at  t*  s  1 
and  x*  -  O.S  as  a  function  of  Pr  number  is  plotted  in 
Fig.  ll .  The  velocity  decreases  with  the  Prandtl 
number.  As  explained  earlier,  this  is  du**  to  the 
decrease  in  temperature  gradient. 


Table  3  Value  of  the  Governing  Parameters  of  the  Five  Cases 


Case 

1 

2 

1 

4 

s 

Re 

12.5 

12.5 

12. b 

12. S 

12.  S 

Pe 

1.25 

0.25 

32. S 

1  .2S 

1  .2S 

Pi 

0.  ’ 

0.02 

0.1 

• .  i 

s 


I ).) 


SS ,000 


5S,lM){i 


10.000 
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more  important.  This  is  reflected  by  the  fact  that  the 
aspect  ratio  (width  divided  by  depth),  as  plotted  in 
fig.  IS,  increases  with  the  Prandtl  number. 

E.  Cooling  Hate 

The  cooling  rate  along  the  interface  of  the  first 
case  is  plotted  in  Fig.  16.  The  cooling  rate  is 
calculated  from  Eq.  (g)  .  The  •  -oling  rate  increases  as 
one  moves  away  from  the  bottom  (center),  then  it 
attends  a  local  maximum  and  decreases,  finally  it 
increases  to  the  edge  of  the  pool.  The  existence  of 
the  local  maximum  is  unusual,  but  expected.  The  reason 
is  because  of  the  existence  of  the  two  counter  rotating 
vortex.  The  cooling  rate  attends  its  local  maximum  at 
the  location  where  the  two  vortex  meet.  In  aU  five 
cases,  the  cooling  rate  at  the  edge  of  the  pool  is 
higher  than  that  at  the  bottom  (center)  of  the  pool. 
Thus,  the  molten  region,  as  it  solidifies,  would  shrink 
faster  in  its  width  than  its  depth.  This  predicted 
trend  is  consistant  with  the  experimental  fact  that  the 
microst ructure  is  finer  at  the  edge  of  the  pool  than  at 
the  bottom  of  the  pool.  The  cooling  rate  at  the  edge 
of  the  pool  as  a  function  of  S  number  and  Pr  number  are 
plotted  m  figs.  1/  and  18,  respectively.  The  cooling 
rate  increases  with  the  S  number  since  higher  the  's 
number,  higher  the  convection  and  thus  more  efficient 
is  the  heat  transfer.  The  cooling  rate  increases  at 
first  and  then  decreases  with  Pr  number.  The  reason 
behind  is  that  there  are  two  parts  m  the  cooling  rate, 
convection  and  conduction.  For  small  Pr  number, 
conduction  is  dominating  so  that  an  increase  m  r*r 
number  decreases  the  conduction  and  thus  the  cooling 
rate.  On  the  other  hand,  where  Pr  number  is  large, 
convection  is  more  important.  Convection  increases 
with  Pr  number,  hence  more  efficient  heat  transfer,  v 


•ij  II  Horizontal  Surface  Velocity 
versus  Prandtl  NuC.ber 


fig.  12  Molten  Region  at.  t*  5  1.0  for  Oifferent  S 
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a  result,  the  cooling  rate  increases. 


Keg  ion  at  t*  1.0  for  Oifferent  pr 
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F.  Solute  Redistribution 

The  recirculating  velocity,  as  predicted  by  this 
model,  is  of  one  or  two  orders  of  magnitude  higher  than 
the  scanning  speed.  This  means  a  fluid  particle  would 
recircu late  several  times  before  it  solidifies.  This 
can  account  for  the  highly  disperse  and  uniform  solute 
redistribution  within  the  molten  region.  The  transport 
process  is  convection  dominated.  The  solute  is 
actually  transported  by  the  carriage  of  the  flowfield, 
molecular  diffusion  play  a  rather  subordinate  role. 
Because  of  the  fact  that  there  are  two  counter  rotating 
vortex,  one  would  expect  to  have  segregation,  if  any, 
at  the  location  where  the  vortex  meet. 

CONCLUSIONS 

This  two-dimensional  convective  heat  transfer  and 
fluid  flow  analysis  has  revealed  many  important  aspects 
of  the  surface  tension  driven  fluid  flow  in  weld  pool 
and  Its  effect  on  pool  shape,  cooling  rate,  velocity 
field  and  solute  redistribution.  The  important 
findings  are  as  follows: 

1.  ftec i rculat f ng  velocity  is  predicted  with  the  flow 
one  or  two  order  of  magnitude  higher  than  that  of 
the  moving  heat  source.  There  are  two  counter 
rotating  vortexes  except  in  the  high  Prandtl 
number  (2.6)  case. 

2.  Surface  temperature  gradient,  which  is  the  driving 
force  for  the  fluid  flow,  is  maximum  at  the  edge 
o*  the  beam  leading  to  the  most  rigorous  outward 

velocity . 

3.  As  Prandtl  number  increases,  the  aspect  ratio 
(width  divided  by  depth)  increases. 

As  surface  tension  number  Increases  from  30,000  to 
hn.OOO,  the  aspect  ratio  increases.  However,  when 
't  m.reases  from  SS.000  to  100,000,  the  aspect 
rit  m  decreases  leading  to  "key-hole''  mode 
rat  1  on  due  to  the  counter  rotating  vortex. 

*>.  ’*«.  (poling  rate  at  the  edge  of  the  pool  is  found 

to  he  higher  than  that  at  the  bottom  of  the  pool 
beHwv  the  centerline.  This  implies  the 
so '  •  1  »  *  i  cat  i  on  will  proceed  from  the  edge  of  the 

po-  » . 


Prandtl  % uriber 


Distonce  Along  interface 


Fig.  16  Cooling  Rate  along  Interfa ■  •>  for  Different  t* 
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uraph  or  cooling  Kate  dt  Pool  tdgt‘ 
ve* sus  Sur  face  Tension  Number 


rig.  18  Graph  or  Cooling  Rate  at  Pool 
Edge  versus  Prandtl  Number 


b.  As  •■’rjrhJtl  number  increases  from  0.02  to  0.1,  tne 
cnjlmg  rate  at  tne  edge  of  the  pool  decreases. 
ii.il  it  increases  *#nen  the  Prandtl  number  increases 
f rc *  0.1  to  2.6. 

7.  The  cooling  rate  at  tne  edge  of  the  pool  increases 
with  the  surface  tension  number. 

8.  Because  of  the  high  reci rculating  flow*  uniform 
solute  redistribution  is  expected.  The  existence 
of  tne  counter  rotating  -vortexes  in  the  pool 
implies  some  segregation,  if  any,  can  be  expected 
at  tne  location  where  the  vortexes  meet. 
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SURFACE  TREATMENT  OF  SPECIAL  STEELS 
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•  .■  ■  -r.t  ■  tee's  each  raving  two  different 

■  1  ‘  •  ■  i.  i •  w*"'o  1  on  a  (.  e  melted  wits  a 

■.■■  ■  •  «  !■.:  variuus  traverse  speeds  ana  beam 

: '  a  ’teters  . 

'«•  experimental  results  are  described,  toge- 
t ter  witti  some  preliminary  mathematical  modelling. 

IMROU.'.' :  )N 

' oe  use  of  a  laser  to  roduce  rapidly  solidi¬ 
fied  structures  in  the  surface  of  different  mate¬ 
rials  has  recently  received  considerable  attention 
This  is  principally  because  it  is  one  of 
the  few  methods  available  for  producing  these  inter¬ 
esting  structures  on  the  surface  of  articles  of 
use  in  engineering.  Such  as  wear  surfaces, 

The  current  level  of  research  is  aimed  at  cate¬ 
gorising  what  structurescan  be  obtained,  with  some 
work  on  detailed  structural  analysis  (3)  and  others 
on  their  engineering  properties  such  as  wear  resis¬ 
tance,  corrosion  resistance,  and  fatigue  resistance 
(<!)■ 

Discussed  here  are  the  structures  obtained  by 
traversing  a  laser  over  four  different  steels. 

EXPERIMENTAL 

A  Control  Laser  2kW  CW  CO?  laser  was  used.  The 
laser  power  output  was  fixed  at  1.7kW  and  passed 
with  various  beam  diameters  over  four  different 
steels  whose  compositions  were:- 


The  steels  were  received  in  two  states,  the  fully 
austenitised  condition  and  the  fully  heat  treated 
condition,  except  for  the  S82  which  was  surface 
carburised  to  around  1mm  depth  and  then  heat  treated 
at  100°C  for  4  hours  or  140°C  for  4  hours. 

The  range  of  traverse  speeds  explored  was  from 
Sivii/s  to  1200mm/s  while  the  beam  diameter  (defocu- 
sed  Gaussian  beam)  was  varied  from  0.2nm  (focus)  to 
4 . 3mm . 

RESULTS  AND  DISCUSSION 

The  results  were  examined  to  establish: 

1)  The  dimensions  of  the  laser  treated  zone 
2}  The  solidification  mechanism 

3)  The  Crystal  structure 

4)  The  extent  of  any  surface  modification  and  upset. 

THE  DIMENSIONS  OF  THE  LOWER  TREATED  ZONE 

The  variation  of  depth  and  width  of  the  melted 
zone  is  shown  in  fig.l  and  2.  The  laser  coupling 
varies  from  deep  penetration  by  "keyhole"  formation 
to  conduction  limited  melting  and  surface  transfor¬ 
mation  hardening  with  no  surface  melting  (5). 

It  is  apparent  from  fig.l  and  2  that  the  prior 
heat  treatment  had  no  effe-t  on  the  melt  zone,  nei¬ 
ther  did  it  change  the  size  of  the  heat  affected  zore. 
There  is  only  slight  variation  between  the  steels 
with  the  S99  being  a  little  more  difficult  to  pene¬ 
trate.  The  hyperbolic  relationship  between  speed  and 


Material 

c 

Si 

Mn 

S 

P 

Ni 

Cr 

V 

Mo 

582 

0.14 

0.23 

0.35 

0.005 

0.005 

4.03 

1.09 

0.26 

0.37 

0.21 

0.68 

0.022 

0.007 

2.46 

0.69 

- 

0.606 

30HM 

0.43 

1.S7 

0.88 

0.006 

0.009 

1  .90 

0.80 

0.70 

0.35 

-pi  3S 

1  .0 

0.29 

0.42 

0.01 

0.007 

0.17 

1 

|  W 

u£> 

0.02 

- 

I  73 


81 


Fiq .  1 : 


Variation  of  deptn  of  ! azer  melted  zone  with 
traverse  speed  for  S135  Steel  at  three  diffe¬ 
rent  heat  treatment  conditions  and  beam 
diameters . 


Theoretical  G.R  ratio  isobars  for  a  weld 
exhibiting  keyholing. 


•T 

Fig.  2: 


WIDTH  0*  IMMHUTCO  ZOWl 


Variation  of  width  of  laser  melted  zone  with 
traverse  speed  for  S62  Steel  at  two  different 
heat  treatment  conditions  and  three  beam 
di ameters . 


Structural  map  of  surface  carburised  58 2 
Steel  for  a  laser  power  of  1720W 

A-  Area  of  possible  retained  austenite 
heat  treatment  1. 

A*  Area  of  possible  retained  austenite 
heat  treatment  2. 

C-  Conduction  region 
C-  Conduction  limited  melt  zone 
T-  Transformation  hardening  with  no  surface 
melting 

N-  No  visible  effects 
P-  Porosity  and  cracking 


9/fi 


Theoretical  r'R  ratio  isobars  for  a  weld 
exhibiting-  molding. 


t rv i  .11  1.  .-.pi  t-i!  by  .imp;.-  i-iiir  .  (  i- 

>  ’  ’  '  rras.inln';  ,  1  :  ■  -  :  -f  V  1;  •  ■  r 

IimJ;  '  tin--  ljiii  I  u  s  i  on  that  ■  >:>n,  Iir.,,  dpui'j  il- 

P/VD  >  60  +  20  J/mrn^ 

cor.ductior  limited  melting  if  60  +  20  p  P/VU  >  2  * 

2 

U/mm  and  transformation  hardening  if?*  1  ;  P/VD 
2 

;  0,7  *  0,4  J/mm  where  P  =  incident  power  W. 

P  beam  diam,  nr>,  y  -  traverse  speed  in  s. 

SOL  I  OIF 1CAT ION  MECHANISM 

Solidification  is  by  a  dendritic,  or  cellular- 
dendritic  growth.  Fig.  3  shows  a  theoretical  cal¬ 
culation  for  a  keyhole  melt  of  the  thermal  gradient 
G  divided  by  the  solidification  rate  R  (G/R  ratio, 
which  has  been  considered  related  to  the  stability 
of  various  solidification  mechanisms  (6).  A  consi¬ 
derable  variation  in  value  is  calculated  particular¬ 
ly  at  the  chill  edge,  giving  the  possibility  of  a 
thin  equiaxed  zone,  leading  to  a  dendritic  core  with 
a  tendency  towards  cellular  dendritic  growth  towards 
the  centre  surface.  Fig.  4  also  for  a  keyhole  melt 
is  a  plot  of  the  cooling  rate  (G.R)  over  the  melt 
zone.  This  parameter  has  been  noted  (6)  to  be  cor¬ 
related  to  the  dendrite  arm  spacing.  This  parti¬ 
cular  calculation  suggests  a  variation  in  arm  spa¬ 
cing  with  a  finer  arm  spacing  at  the 

centre  of  the  keyhole  and  the  central  surface, obser¬ 
ved  values  in  a  comparable  keyhole  run  vary  fromlpm 
toupm.  For  this  calculation  it  was  assummed  that 
the  dendrites  grew  along  the  maximum  thermal  gradi¬ 
ent,  but  examination  of  dendrite  growth  directions 
showed  that  although  basically  following  the  predi¬ 
cted  direction  there  were  notable  exceptions  which 
must  be  due  to  convective  movement  within  the  melt. 

Some  solidification  cracking  and  blind  weld 
porosity  was  observed  in  the  deeper  keyhole  samples. 
Fig.  5  is  a  structural  map  of  the  S82  steel. 

CRYSTAL  STRUCTURE 

Optical  microscopy  showed  martensite  and  a 
white  etching  retained  austenite  .  The  hardness 
figures  were  used  to  confirm  these  structures. 

The  high  C,  low  Si  and  V  steels  favoured  retai¬ 
ned  austenite  at  the  faster  processing  speeds.  The 
300M  steel  did  not  show  any  retained  austenite. 

The  S135  steel  showed  an  unusual  banded  struc¬ 
ture  which  is  apparently  molten  since  the  surface 
had  moved  in  some  samples  but  was  apparently  not 


'h;  i  t » ■ '  .  r  ••  "  r,-  t,Lif,  p,  f , .  T  5  [ji„-  J  :M*-rr  -I  ;•  *  ' 

navv  beei.  partially  molten  due  to  the  large  OiMe 
rentes  between  solidus  and  liquidus  in  this  1  ;C 
steel . 

SURFACE  PROFILE 


By  using  larger  beam  diameters  smooth  suria  ef 
were  obtained. 

CONCLUSION 

data  on  the  dimensions  of  laser  treated  zones 
are  given,  together  with  preliminary  mathematical 
predictions  of  the  solidification  mechanism.  The 
crystal  structure  is  martensite  with  some  retained 
austenite  in  the  high  C  low  Si  and  V  steels. 

REFERENCES 

1.  B.H.  hear,  E.M.  Breinan  L.E.  Greenwald  and 
C.M.  Banas,  Tech,  paper  M.R.  76-867,  1976. 

2.  S.M.  Copley,  D.  Beck,  0.  Esquivel  and  M.  Bass, 
"Laser- Solid  Interactions  and  Laser  Processing" 
pl61,  edited  by  S.D.  Ferris  and  H.J.  Leamy  and 
J.M.  Poate,  American  Institute  of  Physics,  New 
York  1979. 

3.  B.G.  Lewis,  D.  Gilbert,  and  P.R.  Strutt,  "Proce¬ 
ssing  and  Properties  of  High  Speed  Tool  Steels", 
p84,  Conf.  Proc.  AIHE  Las  Vegas  Feb.  1980,  Publ. 
Metal.  Soc.  AIME,  Warrendale  Pa. 

4.  T.R.  Anthony  and  H.E.  Cline,  J.  Appl.  Phys.  p  1 318 
vol .  49,  No. 3,  March  1978. 

5.  G.  Christodoulou  and  W.M.  Steen,  "Laser  ’81  Conf 
Munich  1981  to  be  published. 

6.  M.C.  Flemings  "Solidification  Processing",  p75- 
76  4  146-154  McGraw  Hill,  New  York  1974. 


I  76 


82 


Qpf>e>ol>ix  F 
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ABSTR ACT 

Thr  process  of  laser  cladding  using  gas 
borne  injected  powder  has  already  achieved 
industrial  acceptance.  The  cladding  material 
in  the  form  of  a  fine  powder  is  injected 
into  the  laser  generated  melt  pool  on  a 
moving  substrate.  The  present  model  is 
aimed  at  calculating  the  resulting  clad 
weld  bead  dimensions  and  the  temperature 
distribution  in  the  heated  zone.  A  simula¬ 
tion  model  is  also  presented  for  the  over¬ 
lapping  of  single  clad  tracks  to  produce 
layers  of  uniform  thickness. 


I NTRODUCT I  ON 

Cladding  is  one  of  the  many  material  proces¬ 
sing  applications  for  which  a  laser  has 
been  and  is  being  used.  Fig  1  illustrates 
the  powder  injection  process,  currently 
being  used  by  Rolls  Royce  M). 


fig  1.  Laser  cladding  by  powder  injection. 

The  cladding  material  in  the  form  of  a 
fine  powder  (usually  of  non-spherical  shape 
and  average  si  if  of  70  microns)  1  s  injected 
into  the  laser -generated  melt  pool  on  the 
moving  Substrate  Th  ”  other  r-e  ,sf 


laser  cladding  is  by  melting  a  pre-placed 
powder  layer.  The  use  of  the  cladding 
material  in  the  form  of  a  powder  is  pre¬ 
ferred  as  an  instantaneous  layer  of  powder 
on  or  near  the  surface  appears  to  enhance 
the  laser  energy  coupling  efficiency. 


The  mathematical  model  presented  here 
is  applicable  to  the  injection  process. 
The  main  parameters  of  thi9  process  and 
the  effect  of  these  parameters  on  the 
dependent  parameters  have  been  studied 
and  reported  elsewhere  (2). 


The  model  described  here  assumes  a  quasi¬ 
steady  state  with  the  axis  of  reference 
fixed  relative  to  the  laser  beam,  the 
beam  being  incident  normal  to  the  substrate 
and  stationery  with  respect  to  the  earth. 

The  beam  is  assumed  parallel,  with  a  mathe¬ 
matically  defined  power  intensity  distri¬ 
bution  similar  to  that  produced  by  a 
T.E.M.00  mode.  The  spot  radius  is  defined 
as  the  radial  distance  where  the  central 
intensity  is  reduced  by  a  factor  1 /e . 

The  present  model,  although  not  lacking 
in  generality.  is  aimed  at  co-relating 
experimental  data  obtained  for  flat  plate 
cladding  of  stainless  3)6  to  a  mild  steel 
EN3  substrate,  using  a  2  kw  COi  laser. 

The  temperature  dependence  of  some  physical 
properties  is  allowed  for.  Effects  due 
to  latent  heats  are  accounted.  Heats 

generated  by  chemical  reactions  and  ‘body 
convection*  due  to  melt  pool  turbulence 
are  neglected. 


A  semi -spher ica 1  reflecting  device  was 
used  to  measure  the  total  reflected  p  o  w  e  4- 
trom  a  shot  blasted  mild  steel  surface, 
fig  2.  A  reflectivity  value  of  0.4  was 
obtained  fo  r  a  normal  incident  CO,  1  a  r-  “  r 
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*  beads  of  similar  widths  but  of  different 
hpghts  when  cladded  on  the  same  substrate, 
provided  no  chemical  reactions  occur  (e.g. 
excessive  surface  oxidation  or  formation 
of  inter-metallic  compounds).  Fig  6. 
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iii 


Fig  6.  Similarity  of  substrate  melt  width 
(iii)  and  clad  bead  width  (i,ii). 
Substrate  mild  steel,  clad  -  ( i i ) 
Brass,  (i  )  Stainless  steel. 

Beads  of  cladding  materials  with  a  higher 
melting  point  than  the  substrate,  spread 
out  to  the  full  width  of  the  substrate 
surface  melt.  If  the  energy  input  is  insuf¬ 
ficient  to  generate  temperatures  higher 
than  the  melting  poinc  of  the  clad,  then 
individual  particles  are  observed  entrapped 
in  the  substrate  molten  pool.  Similarity 
of  the  two  widths  will  be  greater  for  clad/ 
substrate  combinations  with  similar  high 
temperature  physical  properties.  Satis¬ 
factory  cladding  is  not  possible  when  the 

melting  point  of  the  clad  is  much  higher 
than  the  substrate  (2). 

(b)  The  section  profile  of  a  single  clad 

bead  is  a  segment  of  a  circle.  This  profile 
has  been  observed  to  be  independent  of 
all  process  parameters.  Clad  beads  of 
stainless  steel,  iron,  brass,  stellite 
and  nickel  were  examined. 

Laser  cladding  is  a  comparat i vel y  calm 

process  relative  to  plasma,  flame  or  arc 

processes.  There  are  no  electro-magnetic 
forces  or  gas  jets  to  disturb  the  molten 

pool.  Also  it  is  evident  that  there  is 

little  molten  pool  turbulence.  This  was 

verified  by  cladding  pure  iron  to  a  stainless 
steel  (18%  Cr)  substrate  under  conditions 
which  would  create  a  high  dilution  level. 
Although  there  was  about  301  dilution, 
(i.e.  30%  of  the  molten  pool  extending 


Fig  7.  Three  basic  clad  bead  section 
files.  Clad  -  Stainless  steel, 
Substrate  -  Mild  steel. 


MODEL  OUTLINE 

Certain  parameters  are  first  calc-slate 
using  an  existing  finite  difference  las.- 
heating  model  (3)  assuming  no  format  if-: 
of  a  clad  bead  and  no  powder  inv*cf.i(in. 
These  parameters  are: 

ic)  Surface  melt  width  and  centreline  melt 
length 

(d>  Surface  melt  area 

(e)  Position  of  the  beam  centre  in  th- 
molten  pool 

(f)  The  length  of  the  heated  rone  in  th*: 
x,y,z  direction:  XR,yR,2R. 

The  objective  of  the  above  calculation 
is  to  obtain  a  first  estimation  of  th* 
bead  section  profile. 

Observation  (a)  allows  a  first  estimation 
of  the  clad  bead  width  to  be  obtained  using 
(c)  above.  A  first  estimation  of  the  clad 
bead  height  is  calculated  by  using  td) 
above,  assuming  that  the  clad  bead  is  made 
from  all  the  particles  falling  on  to  the 
molten  pool  area  and  by  expressing  the 
injected  mass  flow  as  mass  flow  per  unit 
area. 


A  cartesian  matrix  is  set  up  allowing  suf¬ 
ficient  grid  points  inside  the  estimated 
clad  bead  profile.  The  lengths  of  the 
heat  affected  zone  calculated  in  (f)  above, 
are  used  as  initial  limiting  lengths. 


The  power  input  to  each  surface  element 
is  then  calculated  using  a  procedure  similar 
to  that  described  in  (3),  i.e.: 


POWER ( K , I )• 


$■ 


P ( x , y ) dydx 


where  'R*  is  the  region  bounded  by  grid  point 
(K,I)  and  P(x,y)  la  the  power  intensity 
distribution  function. 
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(iv) 


•'  i  g  ft  .  Matrix  set  up.  (  i  )  General  a  r  :  a  -.  g  •  -  • 
(11)  Clad  bead  body  -  trar.  svers" 
section  (iii)  Clad  bead  molten  ^ , 
end  -  longitudinal  section  '  i  v  .  v  : 

'  A  ’  of  (  i  i  i  )  ,  (v)  Macrograph  of  it. 

(vi)  Macrograph  of  (iv). 

The  shadowing  effect  of  the  particle  :luju 
is  evaluated  for  each  surface  grid  point. 

The  progressive  temperature  rise  of  a 
powder  particle  traversing  through  the 
laser  beam  and  falling  on  to  a  particular 
surface  grid  point  is  calculated  for  the 
various  energy  inputs  en  route,  allowing 
for  a  forced  convective  loss  due  to  the 
conveying  gas  flow  and  a  radiative  loss 

In  the  subsequent  heat  balance,  the  heat 
content  of  the  partic  es  is  added  to  the 
respective  surface  grid  points  as  a  body 
convection  term  and  the  power  input  is 
adjusted  according  to  the  shadowing  effect 
of  the  particles. 

If  the  resulting  temperature  distribution 
does  not  satisfy  certain  criteria  which 
depict  the  quasi-steady  state  of  the  molten 
pool,  then  the  clad  bead  width  and  height 
are  varied  until  a  satisfactory  distribution 
is  obtained. 


MATRIX  SET  UP  -  TEMP  (K.r.LJ 

The  circular  bead  section  is  modelled 
as  a  stepped  cartesian  grid  as  shown  in 
Fig  8.  Symmetry  allows  only  one  half 
of  the  bead  to  be  considered.  The  boundary 
of  the  grid  is  defined  by  the  bead  radius. 
The  front  end  of  the  clad  track,  i.e. 
the  molten  pool  end,  is  modelled  as  follows: 
With  reference  to  Fig  8(ili),  the  molten 
pool  elongation  factor  ‘RM*  is  calculated 
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.  s  .1  *  i  .>•:  a*  -  *  ’-  •  .  •*  i  rr  ••  n  t  i  •  s  r  •  !  » '  .  • 
r<  the  molten  p'.-ol.  Th<*  be*-:  ta-J.u:.  : 
calculated  at  each  section  i  e.g.  ^sin-j 
'HI*  as  height  and  2Y1  as  width  at  x  »  X', 
and  a  stepped  grid  synchronous  with  the 
clad  body  grid  is  set  ip. 

Twelve  grid  points  are  allowed  ;  ‘x- 

direction  between  Cl  and  C2 .  T;>#>  total 

grid  size  is  10  v  >  x  ’*5. 


THE  SHADOWING  F  F  F  EOT  it  T  H  K  PARTICLE  CL‘»Ur 
-  FLOSS  ( K ,  I  I 

Assuming  a  weighted  average  particle  dia¬ 
meter,  the  number  of  particles  per  unit 
volume  is  first  calculated.  This  is  a 

function  of  the  particle  velocity,  particle 
diameter,  mass  flow  and  injector  tub*' 
ri  1  are*  e  r  . 

Fig  9  shows  the  powder  cloud  in  relation 
to  the  laser  beam.  The  powder  cloud  is 
considered  to  be  of  square  section  for 
the  purpose  of  statistical  analysis,  a 
correction  being  made  at  a  later  stage. 

The  shadowing  effect  is  evaluated  as  a 
percentage  laser  power  loss  for  each  surface 
element  within  the  powder  cloud.  As  the 

particle  distribution  is  considered  to 
be  uniform,  the  problem  is  reduced  to  a 
single  dimension,  the  effect  being  indepen¬ 
dent  of  the  element  location  in  the  'y1 

direction. 

Reference  to  Pig  9,  the  instantaneous  number 
of  particles  in  the  volume  of  space  above 
each  element  is  calculated.  A  two  dimen¬ 
sional  matrix  is  then  set  up,  the  matiix 
lengths  being  the  element  length  i-.  the 

x-dlrection  and  the  cloud  range  in  the 

y  direction.  The  matrix  positions  are 

spaced  one  particle  diameter  apart.  The 
particles  are  randomly  distributed  into 
the  matrix  postions  and  using  a  statistical 
technique  described  in  M),  the  percentage 
laser  power  loss  is  expressed  as  the  percen¬ 
tage  occupied  space  in  the  matrix  allowing 
for  particles  situated  on  top  of  each  other 
and  a  correction  for  the  area  difference 
between  square  matrix  and  circular  particle 
sect  ion . 


Fig  9.  Powder  clojd  in  relation  t<  o-arr  ar 
substrate  surface. 


Fig  10.  The  shadowing  effect  of  the  particles 
on  the  gaussian  intensity  distributee 
of  the  beam. 


PARTICLE  TEMPERATURES  -  DELT  (K.M 


Pig  10  illustrates  the  shadowing  effect 
of  the  particles  on  the  gaussian  power 
intensity  distribution  of  the  beam.  The 
total  power  loas  is  usually  10-20%  for 
the  range  of  cladding  parameters  used. 


The  energy  inputs  to  a  particle  travelling 
through  the  beam  will  be:  directly  from 
the  laser,  reflected  energy  from  the  sub¬ 
strate  and  reflected  energy  from  other 
particles.  Aa  regards  tha  latter,  the 
multiple  reflections  in  tha  powder  cloud 
are  modelled  as  a  geometric  progression 
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The  traverse  time  across  a  typical  element 
shown  outlined  in  Fig  9,  is  calculated 
assuming  the  particle  travels  along  a 
straight  line  inclined  to  the  substrate 
at  an  angle  equal  to  the  injection  angle. 
The  power  input  to  the  particles  in  the 
vertical  column  containing  thic  element 
will  be  HI.  <  K  ,  »  i  X  f'OWKR  1  K  ,  I  )  .  The  Taction 
absorbed  is  M  .  -  R  F  )  where  R  F  is  the  reduced 

reflectivity  due  to  multiple  reflections. 
The  fraction  absorbed  from  the  beam  reflected 
from  the  clad/substrate  surface  is  also 
calculated  in  a  similar  manner.  The  energy- 
input  to  a  single  particle  in  the  element 
is  calculated  using  the  traverse  time  and 
assuming  that  the  power  input  to  the  vertical 
Column  is  shared  equally  by  all  the  particles 
in  that  column.  Although  the  particles 
at  the  top  are  likely  to  receive  a  higher 
proportion  of  the  laser  energy  directly, 
the  process  is  reversed  for  the  reflected 
beam  from  the  clad/substrate  surface. 

The  net  effect  is  assumed  to  be  an  equal 
power  sharing.  The  temperature  rise  is 
calculated  assuming  the  particle  to  be 
of  uniform  temperature.  A  convective  and 
a  radiative  loss  is  allowed  for.  The  convec¬ 
tive  loss  is  based  on  an  established  empiri¬ 
cal  relationship  (5,6).  A  conveying  gas 
velocity  of  «  m/s,  relative  to  the  particle 
is  used.  Both  losses  are  based  on  an  average 
of  the  temperature  of  the  particle  at  the 
current  position  and  the  immediate  previous 
position.  The  temperature  of  the  current 
position  is  iterated  to  within  10%  of  the 
calculated  current  temperature  uaing  a 
trial  and  error  loop.  The  progressive 

temperature  rise  of  the  particle  is  calcula¬ 
ted  until  arrival  on  the  clad/substrate 
surface.  The  model  calculates  the  highest 
temperature  of  70  micron  stainless  steel 
particles  injected  through  a  5  mm  beam 
of  1030k  to  be  1200k  for  a  powder  mass 
flow  of  0.2  g/a,  particle  velocity  of  1.6  m/s 
and  injection  angle  of  30  deg. 


HEAT  BALANCE 

The  usual  three  dimensional  finite  differ¬ 
ence  heat  balance  equation  is  used  (7,8,9), 
with  two  'body  convection1  terms.  A  body 
convection  term  is  used  for  the  'kinetic 
heat  flow'  in  the  substrate  traverse  direc¬ 
tion,  i.e.  the  'x*  direction.  This  is 
Similar  to  that  described  in  (3).  The 
other  term  accounts  for  the  powder  injection 
proceaa  as  follows.- 


Z_"  ,  >  ^  % 


t  =  temperature  1  - 1  a  part ic:-  , ,  - 

s  1  ze  range  ,  1 . e .  . 

n  *  1  particle  size  SO  ym 

n  =  2  60pm  particle  size  So  Pm 

n  =  k 

m  =  X  M  where  is  the  fr.ii'tioi;  >  ■  f  '  n  •  •  .  -  . 

n  ' 

size  range  in  the  powder  mixture.  ft  j-*j. 
mass  injection  per  unit  time  to  * n  *  7*1  !  ;  • 
surface  area. 

Only  one  size  range  is  cnv;  J  •  , 

we-,  ghted  average  si  ?»  -.  s  .  s  e  d  . 

The  term  is  add*  .1  ,n  l  ,  !  >  t  s  .  :  » 

points  which  are  molten,  •  h  .  .  .  • 

for  other  points  the  part  ;  •  ;  e  s  w .  , 

off*  with  negligible  hta'  ♦  r  a  r.  s  f  *  * 

A  further  assumption  is  that  •  5  »  :  , 

me  Its  and  attains  the  grid  pc,  in'  t  »••*•}•*•  1  ■» »  - 1  »• 
within  the  grid  point  bo-ndar,  F  r  som,- 

larger  particles,  this  may  not  be  true 

as  ejection  of  particles  fr  -  the  mu ’ter 
pool  has  been  observed. 

The  latent  heat  effects  are  jo 
for  using  the  'temperature  e  <j  u  1  v  a  l  •»  r  r  »*  ‘ 

method.  If  a  grid  point  temperature  rea  ;•««**. 
the  melting  point  (l.e.  an  averaqe  be'weer 
solidus  and  liquidus  t empe r a t u r e s  1  then 
the  temperature  is  reduced  by  the  value 


A  convective  and  a  radiative  loss  is  allowed 
for  all  surface  grid  points. 

The  temperature  dependence  of  the  thermal 
conductivity  is  accounted  by  using  linear 
functions  (10,11,12)  for  the  clad  (stainless 
316)  and  the  substrate  (mild  steel  e  N  J  > 
materials,  the  values  being  calculated 
en  route  in  the  iteration  process,  and 
1  imited  at  1650K . 

The  temperature  dependence  of  the  specific 
heat  is  allowed  for  the  mild  steel  sub¬ 
strata.  The  specific  heat  of  stainless 

316  clad  is  assumed  constant  as  variations 
with  temperature  are  not  significant. 

Convergence  is  usually  achieved  after 
500  iteration  cycles.  A  low  relaxation 
parameter  coupled  with  a  low  convergence 
criteria  is  used . 
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(h)  The  temperature  at  point  Cl  li.e.  the 
fully  developed  clad  bead  height)  should 

be  the  melting  point  of  the  clad. 

< i >  The  temperature  of  points  after  C I 
in  the  positive  x  direction  should  be  below 

the  clad  melting  point  if  they  are  within 
the  powder  in  i**ct  ior.  range.  If  this  is 

not  true,  then  the  fully  developed  bead 

height  will  not  oe  at  point  Cl.  The  above 

argument  should  apply  to  all  other  points 
in  the  molten  pool.  However,  on  1  y  (x,o,o) 
points  are  considered  as  they  are  not  affec¬ 
ted  by  the  'step'  in  the  grid. 

(  j  )  A  continuous  interface  zone  at  cne 
melting-point  of  the  substrate  to  facilitate 
a  fusion  bond.  This  is  obtained  if  'g' 
above  is  satisfied. 

(k)  The  substrate  surface  points  adjoining 
the  clad  bead  boundary  (including  the  ellip¬ 
tical  boundary  at  the  front  end)  should 

not  be  molten. 

The  temperature  distribution  is  accepted 
as  a  solution  if  the  above  criteria  are 

satisfied.  In  operating  the  model,  the 
temperature  distribution  of  the  first  esti¬ 
mated  bead  section  profile  is  observed. 
The  temperatures  of  the  critical  points 
of  a  first  estimated  bead  profile  are  as 
f ol lows : 

Width  •  2.5  mm  Height  ■  0.612  mm 
Temperature  at  Cl  -  1  48  7K 
C2  -  1603K 

C3  -  1575K 

The  above  temperatures  are  seen  to  be  too  low 
the  melting  point  being  1650K.  The  limiting 
lengths  of  the  heated  zone,  XR,YR,ZR  are 
now  adjusted  if  required.  The  program 

is  re-run  with  a  new  width  and  height  value 
until  a  satisfactory  temperature  distribution 
is  obtained. 

Temperature  distribution  for  various  widths 
and  heights  were  studied.  It  was  found 
that  only  a  single  combination  of  width 
and  height  satisfy  all  criteria  concurrently, 
with  experience,  it  is  possible  to  arrive 
at  this  solution  with  4-6  runs  of  the  pro¬ 
gram  . 

Temperatures  at  the  critical  points  of 
the  accepted  bead  profile  are  as  follows: 

Width  -  2.36  mm  Height  -  0.669  mm 


(  l  >  Partial  temperature  i  k  )  d  l  s  t  r  i  o  - 1  i 
transverse  section  at  x  «  -  0  .  5  2  6  mm. 


r 

0 . 000 

0-  262 

0.524 

0.  70  7 

1.04?  1  .  3!  1 

i 

0.  000 
0.K9 
3.297 

3 . 4  4  t> 

3-595 

'  650* 

1  6  4^ 

3.743 

151  i 

'  50e 

1  50f- 

'46  9 

1409  *  3?* 

'  -  399 

’  008 

1  004 

994 

976 

95:  9i- 

'5 

*  Temperature  at  0  2 


(  i i  >  Part ial  t  empe r  a  t a i 
transverse  section 


T  3.000 

Z 

0.000 

0.149 

0.2 97  1 856 
0.446  1810 
0.595  1730 
0.743  1628 
1.399  M91 


1855 

1820  1853 

1730  1735 

1627  1625 

1189  '161 


176  3 

1  703  1  649* 

1612  1587 

1166  1143 


1562  1213 

1112  91 3 


*  Temperature  at  C3.  The  width  is  fully 
developed  before  the  height. 

(iii)  Partial  temperature  (k)  distribution  o* 
transverse  section  at  x  *  2.584  mm. 


T 

0.000 

0.262 

0.524 

0.787 

1.049 

1.31 

1  2.2 96  a* 

0.000 

1  650* 

_ 

. 

_ 

. 

0.  149 

1604 

1607 

1548 

- 

. 

. 

. 

0.297 

1547 

1543 

15H 

1471 

. 

. 

. 

0.446 

1487 

1482 

1463 

1441 

- 

. 

. 

0.595 

1425 

1421 

1  408 

1391 

1  381 

- 

- 

0.743 

1  362 

1359 

1349 

1332 

1  307 

1254 

978  Subtt  rat* 

1.399 

1077 

1074 

1065 

1050 

1028 

997 

822 

*  Temperature  at  Cl. 

The  above  solution  is  obtained  for  th* 

f ol lowi ng 

pa  rameter s : 

Laser 

power  1  8  3  0 W , 

Spot  diameter  5  mm,  Powder  flow  0.2  g/s, 
Cladding  speed  6.67  mm/s.  Factor  ’  RM'  1.J, 
Factor  *  RN '  0.26. 

A  clad  bead  of  2.24  mm  width  and  0.622  mr 
height  is  produced  experimentally  for  th< 
above  process  parameters.  The  error  ir 

the  mass  deposition  rate  is  14%. 
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The  basic  technique  of  producing  i  n  1  f  o  r  m 
•-lad  layers  by  laser  cladtlinq  is  by  overlap¬ 
ping  of  single  clad  tracks.  A  me  •  d  e 1  n 
developed  to  simulate  this  process.  The 
model  is  based  on  the  following  experimental 
observations  and  assumptions. 

1  *  •  Th  ••  single  t  r  a  c  *  section  profile  1  r, 
a  segment  of  ^  c  i  r c 1 ° 

•  1  »  The  w i d  t  h  at  a  multiple  ti  •;«  layer 
of  ‘N’  tracks  is  equal  to  w*(N-*tx  where  ‘x‘ 
is  the  transvers*  index  and  ’W’  i  «  the  bead 
width. 

(m>  Equal  a  mo  an*  of  clad  metal  i  .«  deposited 
each  time  a  t  r  a  r x  is  laid. 

Using  the  aoove  simple  ideas,  tent,  rkable 
agreement  with  experimental  data  has  been 
obtained.  This  is  illustrated  in  Fig  11, 
where  the  computed  development  pattern 
of  a  clad  layer  is  compared  with  the  experi¬ 
mental.  The  simulation  technique  is  shown 
m  Fig  >2  and  a  macrograph  of  a  clad 
layer  'leading  edge*  is  shown  in  inset  Fig  1 ) 
The  ’leading  edge*  (i.e.  the  distance  before 
the  clad  thickness  become  uniform)  has 
been  substantially  reduced  by  a  higher 
initial  overlap,  the  model  being  used  to 
calculate  the  required  parameters. 


distribution,  beam  diameter,  surface  ref 
tivity,  values  of  which  can  only  be  *>  v  a  1  u 
to  a  limited  accuracy  using  pre«n>n*  ? 
n  i  q  ■  j  e  s  although  a  d  v  a  r.  *'  e  s  are  h  e  i  n  } 
i  1  )  ,  !  4  ,  I'j  )  .  !n  laser  cladding,  additi 
ambiguous  parameters  are  introduced 

the  powder  injection  system.  if  some  exp 
mental  lata  tr*‘  i  v  a  i  1  a  r*  1  e  ,  then  *  ^  e «;  *'•  p 
is  a:  :  realised  if.) 

an  be  jse.j  r(  predict  effect--  !  ■  : 
o  i  '  ess  parameti'is.  !f  no  Sue'  lit  a 

available.  reasonable  information  -  a  r 
provided  if  high  temperature  data  « :  •  •  « v 

able  for  the  clad/substrate  materials 
if  no  excessive  heats  are  produced  bv  - h 
cal  reactions. 


Mitme*  hm 

Fig  11.  Development  pattern  of  a  clad  layer. 
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Fig  '2.  Technique  used  in  simulation  of  overlapping  single  track  clad  beads. 


END 

DATE 

FILMED 


/" 


DTIC 


1 4 


END 

DATE 

FILMED 


DTIC 


*  f .  ‘  t 


